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ABSTRACT 


I  critically  examine  several  acoustic  radiation  mechanisms 
that  have  been  suggested  to  explain  the  ultra- low- frequency 
emissions  from  certain  severe  thunderstorms  that  are  observed  on 
the  ground  and  in  the  ionosphere.  The  candidates  are:   (a)  simple 
acoustic  sources  (monopole,  dipole)  related  to  the  expansion  of 
heated  storm  air  and  to  the  storm-scale  circulation  pattern;   (b) 
the  random  acoustic  noise  from  "turbulent"  motions  inside  storms; 
(c)  low- frequency  thunder,  i.e.,  the  radiation  from  lightning 
discharges  or  electrostatic  relaxation  inside  storms;   (d)  vortex 
sound,  the  radiation  from  several  kinds  of  unsteady  vorticity  or 
the  interaction  of  multiple-vortex  systems;  and  (e)  thermo- 
mechanical  oscillations  driven  by  the  interaction  between  latent- 
heat  release  and  the  storm's  updraft.  Each  mechanism  is  evaluated 
for  realism  in  terms  of  actual  storm  processes,  and  the  predicted 
acoustic  power  and  its  spectral  content  are  compared  with  those 
required  by  acoustic  measurements.   The  mechanism  most  consistent 
with  the  observations  appears  to  be  a  form  of  vortex  radiation, 
in  which  instabilities  and,  ultimately,  multiple  vortices  form 
about  the  periphery  of  a  larger,  mesoscale  vortex,  and  radiate 
narrow-band  sound  as  they  spin  about  a  common  axis. 


INFRASOUND  FROM  CONVECTIVE  STORMS, 
PART  II:   A  CRITIQUE  OF  SOURCE  CANDIDATES 

T.  M.  Georges 
1.   RATIONALE 


One  of  the  most  persistent  problems  arising  from  observations  of  acous- 
tic-gravity waves  in  the  atmosphere  has  been  identifying  natural  sources  of 
the  waves  and  establishing  their  radiation  mechanisms.   In  many  cases,  empir- 
ical associations  have  been  found  between  certain  kinds  of  waves  and  parti- 
cular geophysical  phenomena,  but  even  so,  the  number  of  completely  satisfac- 
tory models  of  source  mechanisms  is  small  indeed.  Emission  processes  are  so 
elusive  mainly  because  the  kinds  of  direct  atmospheric  observations  we  would 
like  are  not  available  in  sufficient  quantity  or  with  sufficient  resolution; 
consequently,  we  are  usually  forced  to  adopt  rather  indirect  and  inductive 
approaches  in  searching  for  the  correct  mechanisms. 

This  report  is  about  my  search  for  a  satisfactory  model  of  the  infrasound 
that  some  severe  convective  storms  radiate.   I  have  given  necessary  background 
information  in  Part  I  (Georges,  1973),  which  summarized  what  was  then  known 
about  two  different  manifestations  of  the  waves:   one  in  radio  soundings  of 
the  ionospheric  F  region,  the  other  in  recordings  of  microbarometric  pressure 
fluctuations  at  the  earth's  surface.   For  brevity,  that  background  will  not  be 
repeated  here. 

In  discussing  the  plausibility  of  various  possible  mechanisms  with  col- 
leagues, I  have  been  consistently  met  with  responses  like,  "Yes,  but  have  you 
considered  such-and-such  mechanism?",  usually  with  the  implication  that  "their" 
mechanism  is  the  obvious  one  to  consider.   These  encounters  led  me  to  present 
my  analysis  of  the  source  mechanism  as  a  juxtaposition  of  the  many  that  have 
been  suggested.   It  is  not  my  purpose  to  set  up  straw  men  to  be  sacrificed  to 
a  mechanism  that  I  have  known  all  along  to  be  the  correct  one;  when  I  began 
this  analysis,  each  candidate  appeared  to  have  merits  that  could  not  be  dis- 
missed easily,  and  each  of  them  required  quantitative  testing.  The  objective 
of  this  testing  was  to  identify  the  "most  likely"  of  the  proposed  models, 
based  on  a  comparison  with  known  wave  observables  and  some  simple  theoretical 
arguments . 

The  analysis  of  the  several  source  mechanisms  spans  diverse  physical 
disciplines.   This,  in  addition  to  the  requirement  for  a  thorough  understand- 
ing of  severe-storm  dynamics,  suggested  that  this  should  be  a  review  paper 
containing  more  than  the  customary  amount  of  tutorial  material.   At  the  ex- 
pense of  some  redundancy,  I  have  arranged  the  contents  so  that  some  sections 
can  be  skipped  or  passed  over  quickly  by  those  who,  for  example,  don't  need  to 
review  severe- storm  dynamics,  or  who  have  already  decided  that  a  given  mech- 
anism is  not  a  viable  one. 


With  one  exception,  the  processes  I  will  discuss  are  all  aerodynamic, 
that  is,  the  waves  radiate  from  the  free  motions  of  the  air  itself.   It  has 
long  been  suspected  that  the  bulk  of  the  acoustic- gravity  wave  energy  in  the 
atmosphere  arises  from  aerodynamic  sources,  but  the  general  theory  of  aero- 
dynamic sound,  in  its  present  form  (cf.  Lighthill,  1962;  Proudman,  1952; 
Stein,  1967),  is  difficult  to  apply  to  the  specific  atmospheric  processes 
thought  to  be  responsible.   It  is  useful  to  examine  aerodynamic  sound  from 
two  viewpoints:   one  treats  the  emitting  air  flow  stochastically,  that  is,  as 
"turbulence";  the  other  seeks  sound- generating  air  flows  that  can  be  modeled 
deterministically,  or  "exactly."  I  conclude  that  neither  viewpoint  alone 
yields  a  completely  satisfactory  model  and  that  a  hybrid  model  is  required  to 
explain  the  entire  emission  spectrum. 

I  begin  in  the  following  section  by  reviewing  separately  the  particular 
aspects  of  convective- storm  dynamics  that  are  related  to  the  acoustic-emission 
models  I  will  discuss.   Each  of  the  proposed  mechanisms  will  be  introduced, 
but  quantitative  examination  of  each  is  postponed  until  the  succeeding  five 
sections.   Table  1.1  summarizes  the  different  kinds  of  acoustic  sources  to  be 
considered.   Section  8  is  an  assessment  of  the  relative  merits  of  all  the 
models  discussed,  in  which  a  most  reasonable  model  is  suggested.  I  also  pro- 
pose refinements  to  the  model  chosen  and  suggest  some  critical  observational 
tests. 


Table   1.1     Source  Candidates 


1.  Simple  sources  in  cumulus  growth 

A.  Latent  heat  release  (monopole) 

B.  Cell  circulation  (dipole) 

C.  Buoyancy  oscillations  (dipole) 

2.  Turbulence 

A.  Reynolds  stress  viewpoint  (quadrupoles) 

B.  Random  vorticity  viewpoint 

3.  Electrical 

A.  Discharge  shock  (thunder) 

B.  Electrostatic  relaxation 

C.  Pulse  lengthening  (nonlinear  effects) 

4.  Vortex  Sound 

A.  Aeolian  tone  (wake  vortices) 

B.  Vortex  instabilities  (multiple  vortices) 

5.  Heat-driven  Oscillations 

A.  Combustion  oscillations  (periodic  latent  heat 
release) 

B.  Condensation  in  supersaturated  environment 


CONVECTIVE- STORM  DYNAMICS  FROM  A 
WAVE -GENERATION  VIEWPOINT 


The  existence  of  many  reviews  and  monographs  on  storm  dynamics  and  the 
abundance  of  research  papers  they  cite  make  it  unnecessary  to  present  here  a 
self-contained  account  of  the  subject  (see  for  example,  Byers  and  Braham, 
1949;  Newton,  1967;  Atlas  et  al.,  1963).   Instead,  this  section  highlights 
only  those  aspects  of  severe  convective  storms  that  seem  germane  to  the  prob- 
lem of  infrasound  generation:   (a)  the  dynamics  of  the  convection  currents 
inside  storms,  (b)  the  dynamics  of  the  environmental  flow  in  the  storm's 
vicinity,  (c)  the  nature  and  energetics  of  electrical  discharges  in  a  storm, 
and  (d)  energy  release  by  condensation  of  water  vapor.   It  will  emerge  that 
storm  features  with  spatial  scales  of  about  a  kilometer  and  temporal  scales  of 
tens  of  seconds  are  most  likely  responsible  for  the  emissions.  But  because 
details  of  storm  phenomena  on  these  scales  are  virtually  unknown,  direct 
observational  verification  is  not  yet  possible.   Inferences  about  the  smaller 
scale  structure  and  motions  of  storms  are  drawn  mainly  from  what  we  know  about 
larger- scale  features.   The  most  suitable  model  will  therefore  be  judged  on 
the  basis  of  realism,  theoretical  consistency,  and  conformity  with  the  wave 
observations. 


2.1  How  Storms  Begin 

A  convective  or  cumulus  cloud  forms  when  moisture- laden  air  rises  above 
the  level  in  the  atmosphere  where,  because  of  adiabatic  cooling,  it  can  no 
longer  retain  all  of  its  moisture  content  in  vapor  form.   Water  vapor  con- 
denses, usually  on  particulate  (condensation)  nuclei  in  the  atmosphere,  into 
cloud  droplets,  some  of  which  may  continue  to  grow  in  size  and  eventually  fall 
out  of  the  cloud  as  precipitation. 

Meteorologists  have  historically  divided  convective  storms  into  two 
kinds,  based  primarily  on  the  synoptic  or  large-scale  conditions  leading  to 
their  development:   Relatively  small,  isolated  thunderstorms,  sometimes  called 
the  air-mass  type,  (Byers  and  Braham,  1949)  grow  from  the  amalgamation  of  many 
smaller  convective  plumes  initiated  by  solar  heating  of  the  earth's  surface. 
These  storms  involve  only  a  single  air  mass,  as  opposed  to  frontal  storms, 
which  develop  near  the  boundary  between  relatively  cold  and  warm  airmasses. 
Fundamentally,  they  represent  the  release  of  potential  energy,  by  local  over- 
turning, of  a  conditional  instability  in  the  atmosphere.   Conditional  in- 
stability is  created  when  synoptic- scale  processes  cause  heavy  (cool  and  dry) 
air  to  overlie  a  layer  of  lighter  (warm  and  moist)  air. 

Storms  occur  only  when  overturning  gives  rise  to  the  condensation  of  a 
sufficient  quantity  of  atmospheric  water  vapor.   When  condensation  begins,  the 
release  of  the  water's  latent  heat  of  vaporization  drastically  alters  the 
energetics  of  the  overturning;  the  increased  buoyancy  provides  positive  feed- 
back that  amplifies  the  overturning  until  some  external  limit  (such  as  mois- 
ture supply)  stops  it.   When  and  where  air-mass  storms  develop  on  a  synoptic 


scale  are  determined  mainly  by  situations  that  supply  moist  air  at  low  levels, 
but  on  a  local  scale,  their  development  is  controlled  by  nonuniformities  in 
surface  heating,  by  topography  and  by  mesoscale  atmospheric  wave  motions. 
Furthermore,  storm  development  is  inhibited  by  elevated  temperature  inver- 
sions.  Because  these  mesoscale  (10-100  km)  controlling  factors  are  poorly 
resolved  on  the  synoptic  scale,  the  specific  times  and  locations  of  air-mass 
thunderstorm  outbreaks  are  very  difficult  to  predict. 

But  most  very  large  and  severe  thunderstorms  occur  near  synoptic- scale 
fronts  or  boundaries  between  air  masses.   The  classic  pattern  is  typical  of 
the  Great  Plains  in  the  spring  and  summer:   A  tongue  of  warm,  moist  Gulf  air 
pushing  northward  at  low  levels  encounters  a  mass  of  cool,  dry  air  flowing 
south-eastward.   The  giant  Great  Plains  thunderstorms  form  near  the  interface 
(or  front)  between  the  two  air  masses.   Frontal  storms  are  generally  more 
severe  than  the  air-mass  type  because  there  is  additional  dynamic  support 
leading  to  deeper,  more  vigorous,  and  more  organized  convection.   Such  thun- 
derstorms often  form  continually  in  squall  lines  that  parallel  the  front  and 
are  pushed  ahead  of  the  front  as  it  moves.   The  majority  of  storms  that  have 
been  identified  as  strong  infrasound  emitters  appear  to  be  of  this  latter, 
more  intense  frontal  type. 


2.2  How  Storms  Develop  and  Intensify 

The  intensity  or  severity  of  thunderstorms  has  been  found  empirically  to 
depend  on  the  amount  of  moisture  available  at  low  levels,  the  amount  of  con- 
ditional instability  in  the  atmosphere,  and  the  vector  shear  (vertical  gra- 
dient) of  horizontal  winds.   Once  a  dynamical  mechanism  releases  the  insta- 
bility, the  positive  feedback  of  latent-heat  release  during  water  phase  changes 
provides  the  fuel  for  continued  buoyancy  and  growth,  accounting  for  most  of 
the  total  energy  released  in  a  large  thunderstorm.   Following  the  onset  of 
precipitation,  certain  storms  experience  an  almost  explosive  growth  into  the 
great  anvil- topped  cumulonimbus  that  penetrates  the  tropopause. 

Even  at  this  early  stage  in  our  discussion  of  storm  dynamics,  we  can 
envision  two  possible  infrasound- radiation  mechanisms.   Our  experience  with 
the  noise  from  explosions  of  all  kinds  suggests  an  analogous  process,  but  on  a 
longer  time  scale,  associated  with  the  explosive  energy  release  that  causes 
rapid  cloud  growth.   In  other  words,  we  want  to  ask  whether  a  developing 
cumulonimbus  looks  acoustically  like  a  "slow  explosion,"  with  its  radiated 
frequency  and  efficiency  lowered.   To  find  out,  we  can  represent  the  release 
of  latent-heat  energy  in  a  storm  by  an  acoustic  monopole  and  calculate  the 
acoustic  energy  radiated  by  a  local  injection  of  fluid  energy  equivalent  to 
that  released  in  condensation.   This  estimate  will  be  carried  out  in  Section 
3.1.   We  will  conclude  that  the  required  acoustic  power  is  available  from  a 
monopole  mechanism,  but  that  there  is  no  evidence  that  the  mechanism  is  effec- 
tive in  the  required  infrasonic  frequency  range.   Much  of  the  radiation  from 
cumulus  growth  probably  goes  into  atmospheric  gravity  waves. 


The  combination  of  convection  and  the  release  of  latent  heat  in  storms 
also  suggests  another  class  of  sound- emitting  processes  that  could  arise  from 
a  resonant  coupling  between  the  two.   Thermomechanical  or  heat-driven  oscil- 
lations are  well  known  in  certain  combustion  systems,  for  example,  rocket 
engines  and  gas  furnaces.   Rayleigh  (1878)  recognized  that  if  an  oscillating 
heat  source  has  a  component  of  its  rate  of  heat  release  in  phase  with  an 
acoustic  density  oscillation,  then  the  acoustic  oscillation  will  be  amplified. 
For  amplification  to  occur,  the  heat-release  process  must  depend  on  the  pres- 
sure, density,  or  flow  velocity,  and  acoustic  energy  must  be  fed  back  to  the 
heat  source  in  proper  phase.   Feedback  can  be  either  direct  acoustic  feedback, 
as  when,  for  example,  the  heat  source  is  enclosed  in  an  acoustic  resonator,  or 
by  means  of  the  flow  that  passes  through  the  heat  source.   In  Section  7  we 
will  see  whether  any  such  mechanisms  are  relevant  to  infrasound  generation  in 
convective  storms. 


2.3  Air  Circulation  in  Storms 

Even  though,  for  historical  and  forecasting  purposes,  air-mass  and  fron- 
tal may  be  convenient  storm  classes,  they  are  inappropriate  for  the  classifi- 
cation of  specific  kinematic  and  precipitation  structures.   Though  differing 
perhaps  in  size  and  intensity,  air-mass  and  frontal  storms  can  and  often  do  go 
through  similar  stages  of  development  with  nearly  identical  structural  appear- 
ances.  It  is  therefore  more  appropriate  to  use  a  classification  system  based 
on  the  concept  of  a  precipitation  cell.   This  scheme  has  emerged  with  the 
advent  of  weather  radars.   Such  devices  have  the  capability  of  resolving  the 
fine  structure  (y   1  km  scale  sizes)  of  precipitation  content  —  strictly  radar 
reflectivity.   Small  regions  of  relatively  high  reflectivity  are  identified 
with  individual  convective  elements  or  cells  and  constitute  the  basic  building 
block  of  convective  storms.  Two  kinds  of  cells  commonly  referred  to  are  ordi- 
nary cells  and  supercells  (a  term  first  introduced  by  Browning  (1964)  and 
applied  to  severe  storms  consisting  of  a  single  large  cell). 

Most  storms  are  composed  of  one  or  more  ordinary  cells.   These  cells  are 
relatively  short-lived  and  evolve  through  three  stages  identified  by  Byers  and 
Braham  (1949)  as:   the  cumulus  stage  (with  updrafts  alone),  the  mature  stage 
(with  updrafts  and  downdrafts  coexisting) ,  and  the  dissipating  stage  (with 
downdrafts  throughout).   Ordinary  cells  have  a  lifetime  of  about  an  hour,  with 
a  mature  stage  of  only  10  to  15  minutes.  Thunderstorm  complexes  can  be  com- 
posed of  several  such  cells  at  various  stages  of  development  (Chalon  et  al., 
1976).   A  disorganized,  nearly  random  collection  of  these  cells  typifies  the 
air  mass  thunderstorm.   An  organized  succession  of  cells,  on  the  other  hand, 
is  more  typical  of  frontal  storms.   A  recent  kinematic  model  of  an  organized, 
multicell  hailstorm  has  been  presented  by  Browning  et  al .  (1976)  and  is 
shown  in  figure  1. 
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Figure  1.      Schematic  model  of  an  organized  multicell  hailstorm,   showing 
a  vertical  section  along  the  storm's  direction  of  travel  through  a  se- 
quence of  evolving  cells.      Solid  lines  are  streamlines  of  flow  rela- 
tive to  the  moving  system;    they  are  shown  broken  on  the   left  side  of 
the  figure  to  represent  flow  into  and  out  of  the  plane  and  on  the  right 
side  of  the  figure  to  represent  flow  in  a  plane  a  few  kilometers  closer 
to  the  reader.      The  cells  are  numbered  as  follows:     n+1  =  growing 
cumulus;   n  =  first  radar  echo;   n-1  =  mature;   n-2  =  decaying   (after 
Browning  et  al.3    1976). 


In  the  presence  of  large  conditional  instability  in  the  surface  layer 
capped  by  relatively  dry,  potentially  cold  air  in  the  middle  levels  and  ac- 
companied by  strong  wind  shear  (usually  veering  with  height  of  the  horizontal 
winds),  a  much  more  intense  and  longer  lived  convective  element,  the  super- 
cell,  will  develop.   The  important  features  of  a  supercell  which  travels  to 
the  right  of  the  middle  tropospheric  winds  are  described  by  Browning  (1964) 
and  summarized  in  Figure  2.   Such  a  storm  often  exists  in  a  quasi-steady  state 
for  one  to  several  hours. 


In  view  of  the  frequency  range  of  the  acoustic  waves  whose  source  we  seek 

(tens  to  hundreds  of  seconds  wave  periods)  we  now  focus  our  attention  on  the 

shorter  term  fluctuations  in  storm  airflow  instead  of  on  the  longer  scales 
appropriate  to  a  storm's  life  cycle. 
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Figure   2.      Plan  view   (upper)  and  vertical  section   (lower)   of  a  model  of  a 
severe   local  storm  that  travels  to  the  right  of  middle-level  tropospheric 
winds.      Stippling  in  the  plan  view  is  proportional  to  precipitation  in- 
tensity;  arrows  indicate  the  main  features  of  the  airflow  believed  to  be 
common  to  all  such  storms.      The  circled  T  shows  the  most  likely  place 
for  tornado  formation.      (After  Browning^    1964). 


There  is  substantial  evidence  that  cumulus  growth  in  convective  storms  is 
not  always  steady  but  sometimes  occurs  in  bursts  (cf.  Letzman,  1930;  Workman 
and  Reynolds,  1949;  Abe,  1923;  Priestley,  1953).   In  multicellular  storms,  the 
successive  emergence  of  individual  convective  cells  at  cloud  top  in  a  succes- 
sion of  visual  cloud  towers  often  presents  a  pulsating  appearance  (Scorer  and 
Ludlam,  1953) .   Anderson  (1960) ,  in  a  photographic  study  of  five  orographic 
cumuli  in  Arizona,  found  clear  evidence  of  a  10-min  periodicity  in  cell  growth, 
superimposed  on  which  were  smaller,  shorter  period  "turbulent"  fluctuations. 
He  concluded  that  cumulus  growth  is  a  superposition  of  an  oscillatory  type  of 
cellular  motion  upon  a  steady  growth,  even  for  smaller  cumuli  that  do  not 
develop  into  thunderstorms. 

Newton  (1966)  has  used  radar  to  study  the  pulsations  of  the  cumulus 
towers  of  an  Oklahoma  squall- line  storm  that  penetrated  the  stratosphere.  He 
explains  the  observed  alternate  inflation  and  deflation  of  towers  as  the 
damped  oscillation  of  an  updraft  parcel  as  it  experiences  negative  buoyancy  as 
it  overshoots  into  the  stable  stratosphere  and"then  a  positive  buoyancy  as  it 


returns  to  the  troposphere,  having  been  warmed  by  mixing  with  environmental 
air.   The  oscillation  is  damped  by  form  drag  and  entrainment  of  environmental 
air.   Judging  from  his  sequence  of  radar  storm  profiles,  the  oscillation 
period  is  roughly  15  min. 

Turner  (1966)  has  suggested  that  buoyancy  oscillations  are  controlled 
more  by  evaporative  cooling  at  cloud  top,  which  reverses  the  buoyancy  of  a 
thin  cloud- cap  layer.   His  laboratory  experiments  with  buoyant  plumes  of  salt 
solution  reveal  violent  and  regular  oscillations  as  mixing  with  their  environ- 
ment reverses  their  buoyancy.   For  typical  cloud  parameters,  he  predicts  an 
oscillation  period  of  15  min. 

Note  that  two  different  kinds  of  pulsating  motion  in  cumulus  growth  are 
now  proposed:   One  involves  oscillations  about  the  level  of  neutral  thermal 
buoyancy  of  individual  air  parcels  (which  probably  experience  mixing  with  the 
environment,  however);  the  other  involves  the  successive  (non-oscillatory) 
emergence  of  separate  larger  convective  cells  at  cloud  top.   The  former  (oscil- 
latory) mode  should  occur  to  some  extent  within  every  buoyant  element,  and 
thus  constitutes  a  small-scale  (turbulent)  motion  superimposed  on  nearly  all 
convective  storms.   The  latter  mode,  however,  is  characteristic  of  non- steady 
multicellular  storms. 

The  distinction  between  steady  and  non-steady  storms  is  not  precise;  in 
fact,  at  different  times  in  the  life  of  a  storm,  either  steady  or  pulsating 
modes  of  convection  may  occur.   Newton  (1963)  says  that  a  pulsating  convection 
mode  appears  to  be  favored  in  growing  cumuli  before  downdrafts  develop,  where- 
as forced  convection,  in  which  strong  downdrafts  play  a  lifting  role,  appears 
to  favor  steady  flow.   A  significant  number  of  storms  may  be  of  an  inter- 
mediate character,  that  is,  having  a  "pulsating  steady"  updraft  in  the  sense 
that  a  connected  column  of  rising  air  has  bursts  of  more  intense  convection 
superimposed  on  it.   This  sort  of  behavior  seems  to  be  the  basis  of  the  sug- 
gestion by  Pierce  and  Coroniti  (1966)  that  buoyancy  oscillations  in  growing 
cumuli  act  like  a  piston  to  generate  acoustic- gravity  waves  that  reach  the 
ionosphere  with  detectable  amplitude. 

All  of  these  facets  of  a  time- varying  circulation  pattern  in  convective 
storms  suggest  a  single  acoustic  radiation  model:   an  acoustic  dipole.  An 
incompressible  toroidal  flow  resembling  th'e  central  updraft  and  surrounding 
downdraft  of  a  single  convection  cell  can  be  reproduced  by  the  velocity  field 
of  a  fluid-dipole  distribution  on  a  horizontal,  disk- shaped  surface,  with 
fluid  sources  on  top  and  sinks  on  the  bottom.   Changes  in  the  strength  of  such 
a  dipole  source  would  be  communicated  instantly  to  all  parts  of  an  incompres- 
sible medium.   In  a  compressible  flow  with  low  Mach  number,  the  velocity  field 
of  such  a  source  remains  virtually  the  same  as  in  incompressible  flow,  but 
changes  in  source  strength  radiate  into  the  medium  as  acoustic  energy.  Thus, 
by  representing  fluctuations  in  updraft  strength  in  convective  storms  as  the 
field  of  a  distribution  of  varying  acoustic  dipole  sources,  we  could  estimate 
the  radiated  acoustic  energy  and  its  spectral  content.  The  same  acoustic 
problem  can  be  alternatively  formulated  in  terms  of  changes  in  the  strength  of 
a  ring  vortex,  or  as  the  radiation  from  a  vertically  accelerating  air-mass. 
In  our  calculations,  we  use  the  last  method. 
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The  oscillations  we  have  been  discussing  need  not  be  monochromatic.  It 
has  been  suggested  that  an  air  parcel  undergoing  buoyancy  oscillations  at  the 
local  Brunt  frequency  might  emit  compression  waves  at  twice  that  frequency 
because  the  oscillating  parcel  displaces  adjacent  air  back  and  forth  hori- 
zontally twice  for  every  vertical  parcel  oscillation.   Furthermore,  vertical 
variability  of  the  local  Brunt  frequency  could  cause  non-sinusoidal  oscilla- 
tions with  high  harmonic  content. 


2.4  The  Role  Environmental  Winds  Play 

Early  concepts  of  thunderstorm  development  held  that  strong  vertical  wind 
shear  diminishes  a  storm's  intensity  by  destroying  its  updraft.  However, 
following  the  Thunderstorm  Project  in  1948  and  1949,  researchers  soon  estab- 
lished a  strong  correlation  between  storm  severity  and  strong  middle- level 
winds.   Specifically,  Fawbush  et  al.  (1951),  after  studying  many  tornado- 
spawning  storms  and  their  environments,  found  as  a  necessary  condition  on  the 
synoptic  scale  the  presence  in  the  middle  troposphere  of  a  jet  of  strong  winds 
exceeding  35  knots  (18  m  s-l) .   Strong  middle-level  winds  also  consistently 
accompany  the  most  severe  and  long-lived  hailstorms.   Newton  and  Newton  (1959) 
established  the  same  behavior  for  intense  thundery  rainstorms.   Ludlam  (1963, 
his  Fig.  5)  shows  a  remarkable  coincidence  between  the  geographic  distribution 
of  severe  hailstorms  and  thundersqualls  and  the  mean  position  of  the  500  mb 
jet  axis  during  the  months  in  which  the  storms  occur.   On  the  basis  of  this 
evidence,  the  existence  of  a  strong  middle- level  jet  is  now  regarded  as  an 
essential  ingredient  for  transforming  a  short-lived  thunderstorm  to  one  that 
persists  for  up  to  12  hours  or  more  and  may  travel  several  hundred  miles. 
Perhaps  such  a  jet  is  the  ingredient  that  also  transforms  an  ordinary  thunder- 
storm into  an  infrasound  emitter. 

The  reason  given  for  storm  intensification  by  wind  shear  was  that  a  storm 
exposed  to  changing  winds  at  different  levels  has  its  supply  of  moist  air 
continually  replenished,  whereas  a  storm  in  a  stagnant  environment  soon  runs 
out  of  "fuel". 

In  an  effort  to  quantify  the  dynamical  interaction  between  a  storm  and 
its  environment,  Newton  and  Newton  (1959)  constructed  a  cumulonimbus  model  in 
which,  for  the  first  time,  wind  shear  plays  the  role  of  intensifying  and 
prolonging  convection,  and  causes  continual  cell  regeneration  on  the  storm's 
forward  right  flank.   Because  strong  convection  effectively  destroys  vertical 
shear  within  the  cloud  (and  not  vice  versa) ,  the  incloud  air  assumes  a  hori- 
zontal velocity  intermediate  between  the  ambient  wind  at  upper  and  lower 
levels.   Through  its  ability  to  resist  environmental  shear,  the  storm  now  acts 
as  an  obstacle  in  the  flow  and  dynamic  pressures  are  exerted  on  the  storm 
boundary,  positive  on  the  upstream  side  and  negative  on  the  flanks  and  downwind 
side  (Fujita  and  Arnold,  1963).   These  dynamic  pressures  convert  a  portion  of 
the  mean  flow's  kinetic  energy  to  kinetic  energy  of  the  storm's  vertical 
motions  as  discussed  in  detail  by  Newton  (1967) . 


A  possible  relevance  of  the  storm- environment  interaction  to  sound  genera- 
tion lies  in  the  role  of  large  storms  as  obstacles  in  the  environmental  flow. 
An  analogy  with  our  everyday  experience  is  again  suggested:  wind  blowing 
around  obstacles  such  as  buildings,  trees,  and  wires  produces  familiar  sounds. 
The  nature  of  these  sounds  is  determined  by  the  structure  of  wakes  generated 
in  the  flow  by  the  obstacles.   We  want  to  know  whether  severe  storms  also 
produce  wakes  and,  if  they  do,  whether  low- frequency  sound  is  emitted. 

Fujita  and  Arnold  (1963)  have  reported  wind- field  measurements  in  the 
vicinity  of  a  Kansas  thunderstorm  that  clearly  show  environmental  wind  maxima 
on  the  flanks  of  the  strong-echo  core  of  the  storm,  in  direct  analogy  with  the 
flow  pattern  observed  in  laboratory  flows  past  cylindrical  obstacles.   As 
early  as  1962,  Fulks  hypothesized  the  existence  of  thunderstorm  wake  vortices 
and  a  possible  relation  to  tornado  generation,  but  little  evidence  has  emerged 
to  support  an  association  with  tornadoes.  Even  if  wake  eddies  are  not  directly 
involved  in  tornado  genesis,  there  is  evidence  that  significant  amounts  of  the 
momentum  of  the  mean  flow  are  converted  into  local  concentrations  of  vorticity 
in  or  around  a  storm,  and  at  least  contribute  one  ingredient  favoring  tornado 
development.   Recent  radar  observations  (Jessup,  1972;  Brown  and  Crawford, 
1972;  Lemon,  1976)  confirm  the  role  of  a  thunderstorm  as  an  obstacle  in  the 
flow  and  suggest  that  the  wake  of  a  storm  updraft  contains  vortex  tubes  of 
cumulus  scale.   But  none  of  these  observations  tells  us  much  about  the  dyna- 
mics, the  life  cycles  or  the  ultimate  fate  of  wake  vortices  or  other  details 
of  cumulus  wakes  that  may  be  invisible  to  radars  because  they  exist  outside 
the  storm  cloud. 

One  possibility  is  that  some  vortices  break  free  of  the  storm-obstacle, 
travel  downstream  like  ordinary  wake  eddies,  and  decay  according  to  ordinary 
two-dimensional  models  for  the  diffusion  of  vorticity  (cf.  Wilkins,  1970). 
Alternatively,  the  same  instability  that  causes  the  storm  itself  may  play  a 
role  in  sustaining  convection  in  shed  eddies,  possibly  causing  them  to  develop 
into  self-sufficient  storms  (cf.  Fujita  and  Grandoso,  1968).   A  third  possi- 
bility is  that  wake  eddies  remain  bound  to  the  parent  storm  and  become  inte- 
gral parts  of  its  circulation,  as  suggested,  for  example,  by  Jessup  (1972)  and 
by  Eagleman  (1971).   In  this  third  case,  it  may  be  less  useful  to  view  the 
phenomenon  as  a  wake  eddy  than  simply  as  a  portion  of  the  environmental  flow 
that  becomes  part  of  a  helical  downdraft,  perhaps  intertwining  with  an  up- 
draft, as  sketched  in  the  paper  by  Browning  (1964) .   Lemon  (1976)  ,  for  ex- 
ample, has  suggested  that  some  wake  vortices  may  be  related  to  the  "starting 
vortices"  of  airfoil  theory,  balancing  the  equal  but  opposite  vorticity  in  a 
storm's  rotating  updraft. 

Consequently,  we  are  not  sure  whether  to  model  storm  wakes  as  fully 
turbulent  or  as  regular  Karman  vortex  trains.   Both  kinds  of  wakes  can  emit 
aerodynamic  sound,  but  different  approaches  are  used  to  analyze  them.   Periodic 
wakes  can  best  be  described  in  terms  of  a  regular  vorticity  distribution, 
whereas  fully  turbulent  wakes  have  to  be  modeled  stochastically. 

A  special  case  of  wake-generated  sound  is  the  aeolian  tone,  in  which  the 
Reynolds  number  of  the  flow  is  in  the  proper  range  for  the  wake  to  be  composed 
of  a  highly  periodic  train  of  vortices  that  have  been  shed  from  the  obstacle. 
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Nearly  monochromatic  sound  is  emitted  at  the  eddy-shedding  frequency.   The 
aeolian  harp  and  the  singing  telephone  wire  are  classic  examples.   That  sim- 
ilar vortex  trains  exist  on  mesoscale  and  larger  has  been  confirmed  by  satel- 
lite photographs  of  the  wakes  behind  certain  islands  (Chopra  and  Hubert, 
1964a,  1964b;  Zimmerman,  1969;  Chopra,  1973).  Either  the  warm-air  bubble  that 
sometimes  forms  over  an  island  or  its  mountains  diverts  the  flow  horizontally 
around  the  island  when  a  suitable  temperature  inversion  prevents  much  vertical 
flow  diversion.   The  same  inversion  often  provides  a  cloud  layer  that  makes 
the  resulting  flow  perturbations  visible  on  satellite  photographs.   Similarly, 
as  a  storm  approaches  the  tropopause,  diversion  of  the  environmental  flow 
should  become  mainly  horizontal  because  vertical  deflection  is  inhibited  by 
the  static  stability  of  the  stratosphere. 

In  Section  6.3.2  I  estimate  the  sound  that  could  be  radiated  by  wakes 
that  thunderstorms  might  produce.   I  conclude  that  the  aeolian-tone  mechanism 
could  radiate  the  required  acoustic  power,  but  the  sizes  of  the  wake  eddies 
that  seem  to  exist  would  radiate  at  frequencies  lower  than  those  observed. 


2.5  Vorticity  Inside  Storms 

The  most  significant  single  clue  to  the  origin  of  severe-weather  infra- 
sound is  its  close  association  with  tornadic  storms  in  the  Midwestern  U.S. 
Specifically,  in  a  recent  study  of  infrasound  received  at  Boulder,  Colorado, 
during  an  entire  storm  season,  45  of  the  86  storms  that  were  identified  as 
infrasound  sources  produced  documented  tornadoes  or  funnel  clouds  (Georges, 
1973;  Georges  and  Greene,  1975) .   The  45  tornadic  storms  were  a  small  fraction 
of  all  the  storms  that  occurred  in  the  Midwest  during  the  observing  period, 
yet  they  constituted  a  majority  of  the  infrasound  emitters.  Furthermore,  it  is 
easy  to  imagine  various  difficulties  with  the  tornado-reporting  system,  such 
that  tornadic  (or  at  least  pre- tornadic)  activity  was  indeed  present  in  many 
or  even  all  of  the  remaining  41  infrasound- emitting  storms,  but  was  not  ob- 
served or  reported  for  any  of  several  reasons  (unpopulated  area,  nighttime, 
tornado  obscured  by  cloud  or  precipitation) .  Because  of  this  compelling  ob- 
servational evidence,  storm  mechanisms  associated  with  tornado  production,  and 
specifically  those  that  give  rise  to  concentrations  of  vorticity  in  storms, 
deserve  special  attention  as  potential  infrasound  generators. 

There  is  no  a  priori  reason  for  believing  that  only  fully  developed 
tornadoes  emit  infrasound.   Indeed,  comparisons  of  the  details  of  timing  and 
duration  of  observed  infrasound  with  tornado  reports  suggest  that  such  is  not 
the  case.   Although  observations  are  sparse,  there  is  evidence  for  the  exis- 
tence of  a  larger  class  of  storm- associated  vortices,  some  of  which  are  re- 
lated to  tornado  genesis,  but  others  of  which  grow  and  die  unseen,  either 
entirely  within  a  storm  cloud  or  in  the  clear  air  surrounding  it.   It  is  this 
broader  class  of  vortex  phenomena  whose  dynamics  and  infrasound- generating 
potential  we  now  examine. 
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Portions  of  large  convective  storms  exhibit  rotation  great  enough  to  be 
observed  visually  (especially  in  time- lapse  photography)  and  in  their  radar 
echoes.   The  nature  of  this  kind  of  rotation  is  evident  from  the  flow  pattern 
depicted  in  Figure  2;  the  low  level  inflow  makes  a  roughly  270°  turn  as  it 
convects  upward  through  the  storm  and  is  expelled  at  the  top.   Brown  et  al . 
(1971)  have  suggested  that  the  "hook"  echoes  seen  on  weather  radars  are  prob- 
ably distortions  of  the  precipitation  echo  by  such  rotating  updrafts.   The 
relevance  of  hook  echoes  to  tornado  genesis  is  thus  reasonable  (Fujita,  1973). 
Rotation  in  storm  downdrafts  has  also  been  suggested  (Eagleman,  1971),  but  the 
evidence  for  its  existence  is  less  substantial. 

Doppler  radars  have  revealed  the  existence  of  mesoscale  vortices  (whose 
solid  rotating  cores  are  5  km  or  so  diameter)  inside  some  tornadic  storms 
(Kraus,  1973;  Brown  et  al.,  1973)  as  well  as  in  non-tornadic  storms  (Figure  8 
and  Kropfli  and  Miller,  1975,  1976) .   Evidently,  such  rotations  are  character- 
istic of  many  kinds  of  convective  storms. 

As  this  and  the  preceding  section  suggest,  there  is  no  shortage  of  indi- 
cations (both  theoretical  and  observed)  that  storm-scale  vorticity  exists  in 
profusion  inside  and  in  the  vicinity  of  thunderstorms.   Yet,  on  the  next  scale 
smaller  than  that  of  the  storm  itself,  that  is,  less  than  a  kilometer,  prac- 
tically no  information  about  vorticity  distribution  exists.  These  are  the 
scales  of  motion  that  are  the  most  likely  sources  of  the  observed  acoustic 
radiation. 

We  have  two  choices  with  regard  to  modeling  these  motions  and  estimating 
their  acoustic  radiation:   First,  in  the  absence  of  a  deterministic  spatial 
and  temporal  description  of  the  vorticity  field,  we  could  resort  to  a  sto- 
chastic description  and  use  measured  values  of  the  turbulent  velocity  fluctua- 
tions inside  storms.   The  acoustic  radiation  from  homogeneous,  isotropic 
turbulence  has  been  modeled,  for  example,  by  Lighthill  (1952,  1954,  1962), 
Proudman  (1952)  and  Morse  and  Ingard  (1968).   This  method  only  marginally 
meets  the  radiated-power  requirements  imposed  by  the  observations;  further- 
more, some  of  the  assumptions  about  spatial  and  temporal  stationarity  implicit 
in  the  turbulent  model  are  not  satisfied  by  storm-related  air  motions. 

Our  second  choice  is  to  calculate  the  acoustic  radiation  from  simple 
time-varying  vortex  configurations.   As  examples,  I  chose  a  corotating  vortex 
pair  and  the  vortex- shedding  process  responsible  for  aeolian  tones.  The  aeo- 
lian-tone  mechanism  could  produce  the  required  radiated  power,  but  again,  the 
frequency  range,  for  realistic  storm  parameters,  is  too  low.  I  find  that  the 
corotating  vortex  mechanism  could  be  effective  at  the  observed  frequencies  and 
could  emit  the  required  power  under  conditions  observed  in  the  most  severe 
storms.   Furthermore,  by  invoking  a  multiplicity  of  such  sources  in  various 
stages  of  growth  and  decay,  the  nonstationary  nature  of  the  observed  radiation 
could  be  duplicated. 
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2.6  Thunderstorm  Electricity 

The  exact  mechanism  of  charge  separation  and  the  details  of  electric 
charge  distribution  in  thunderstorms  are  still  under  active  debate.   The  most 
commonly  suggested  mechanisms  suppose  that  charges  of  opposite  sign  attach  to 
particles  that,  because  of  size  differences,  move  at  different  velocities 
relative  to  the  rising  air  stream  of  a  storm's  updraft.   When  a  large  enough 
potential  difference  is  built  up  between  different  parts  of  a  cloud,  or  be- 
tween a  cloud  and  the  ground,  an  electrical  breakdown  of  the  intervening  air 
takes  place  in  a  series  of  steps  that  culminates  in  a  massive  electrical 
discharge  arc. 

Fortunately,  the  aspects  of  thunderstorm  electricity  that  could  be  rele- 
vant to  infrasound  generation  do  not  seem  to  depend  critically  on  these  de- 
tails.  The  reader  interested  in  the  subject  should  consult  the  books  by  Uman 
(1969)  and  Chalmers  (1967)  and  the  many  references  they  cite.  We  will  confine 
our  attention  to  events  associated  with  the  electrical  (lightning)  discharge. 

The  obvious  relevance  of  lightning  to  infrasound  generation  lies  in  the 
possibility  that  the  infrasound  is  simply  low- frequency  thunder.   But  the 
results  of  numerous  measurements  of  the  acoustic  spectrum  of  thunder  (Few 
et  al.,  1967;  Bhartendu,  1971;  Few,  1969;  Holmes  et  al.,  1971)  indicate  that 
there  is  not  nearly  enough  acoustic  energy  in  the  infrasonic  portion  of  the 
spectrum  to  explain  the  observed  infrasound  intensities.  Only  the  possibility 
that  the  spectral  content  of  thunder  may  be  shifted  into  the  infrasonic  range 
by  nonlinear  pulse-stretching  in  the  upper  atmosphere  prevents  us  from  imme- 
diately discarding  thunder  as  a  reasonable  source  candidate.   In  Section  3.3  I 
examine  the  pulse- stretching  mechanism,  and  the  evidence  that  it  actually 
occurs  in  the  atmosphere,  to  see  if  it  could  quantitatively  explain  the  ob- 
served properties  of  severe-weather  infrasound. 

I  also  examine  a  suggestion  by  Dessler  (1973)  that  electromechanical 
relaxation  inside  thunderstorms  following  lightning  discharges  radiates  de- 
tectable infrasound.   The  basic  idea  is  that  in  charged  regions  of  a  thunder- 
cloud, charged  water  droplets  repel  each  other  and  set  up  a  region  of  locally 
reduced  pressure.   When  the  charge  distribution  is  suddenly  neutralized  by  a 
lightning  discharge,  a  sort  of  implosion  occurs,  radiating  a  rarefaction  wave 
whose  dominant  frequency  and  intensity  are  determined  by  the  size  and  shape  of 
the  emitting  region.   I  raise  several  questions  about  the  validity  of  such  a 
model,  and  make  more  accurate  estimates  of  the  radiated  acoustic  energy.   I 
conclude  that  the  relaxation  process,  if  it  occurs,  probably  radiates  infra- 
sonic pulses,  but  that  it  is  probably  not  the  source  of  the  severe-storm 
infrasound. 

The  following  sections  examine  in  detail  the  acoustic  source  mechanisms 
introduced  in  this  section. 
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3.   SIMPLE  ACOUSTIC  SOURCES 

An  acoustic  monopole  is  used  to  model  the  acoustic  radiation  caused  by 
local  injections  of  fluid  mass  or  energy.   An  acoustic  dipole  is  used  to  model 
the  acoustic  radiation  from  local  changes  in  fluid  momentum.   An  acoustic 
quadrupole  is  used  to  model  the  acoustic  radiation  from  local  changes  in  fluid 
momentum  flux.   I  now  examine  the  appropriateness  of  each  of  these  simple 
acoustic  sources  for  describing  processes  known  to  occur  within  thunderstorms. 


3.1  Monopole  Source 

Unless  a  fluid  or  energy  source  actually  exists  in  a  fluid,  monopole 
radiation  is  very  weak  and  any  effective  monopole  radiation  present  is  due  to 
retardation  effects  in  higher  order  sources  that  are  not  small  compared  with  a 
wavelength  (Morse  and  Ingard,  1968).   In  convective  storms,  energy  is  released 
inside  the  cloud  by  water- vapor  condensation  and  sublimation.   At  the  acoustic 
wavelengths  of  interest  here  (tens  of  kilometers) ,  this  source  must  account 
for  virtually  all  monopole  radiation  from  storms.   We  can  make  crude  estimates 
of  the  total  energy  released  and  its  duration,  but  spatial  and  temporal  de- 
tails of  the  release  are  not  known.  Therefore,  we  can  estimate  the  average 
rate  of  acoustic  energy  radiation  but  cannot  yet  know  anything  at  all  about 
the  higher- frequency  (shorter  wave  period  than  the  storm  lifetime)  radiation 
that  may  be  caused  by  fluctuations  in  the  energy- release  process. 

A  crude  but  sufficiently  accurate  lower  bound  on  the  latent  heat  released 
in  a  severe  rainstorm  can  be  obtained  by  calculating  the  heat  of  vaporization 
of  a  storm's  total  rainfall.   Suppose  an  average  5  cm  of  rainfall  from  a 
severe  thunderstorm  covers  a  ground  area  10  x  10  km.  The  total  water  volume 
is  5  x  10"  m^,  and  its  mass  is  5  x  10^  kg.   Condensation  releases  2.5  x  10^ 
J  kg-1,  so  the  total  heat  released  is  1.3  x  lO1^  J. 

This  heat  energy  dh  (per  unit  mass)  is  partitioned  into  (a)  increasing 
the  storm's  internal  energy  (temperature)  and  (b)  doing  work,  mainly  expanding 
the  air  inside  the  storm  (although  small  fractions  go  into  changes  in  its 
kinetic  and  potential  energies) .  Therefore 

dh  -  C  dT  +  pdv,  (3.1) 

where  v  is  the  specific  volume  of  the  gas,  and  C  is  its  specific  heat  at 
constant  volume.  If  we  assume  that  this  process  takes  place  under  nearly 
constant  pressure,  then  we  also  have 

dh  *  C  dT.  (3.2) 

p  *•   > 

We  can  combine  (3.1)  and  (3.2)  to  find  out  what  fraction  of  dh  goes  into 
expanding  the  air: 
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c 

pdv  =  dh  (1  -  ^  ).  (3.3) 

P 

Recognizing  that  C^/C     =   1/y  and  letting  y  =  1.4,  we  obtain,  for  the  change  of 
storm  volume  dV,     " 

pdV  =  0.29dH.  (3.4) 

(Lower  case  letters  refer  to  specific  quantities  per  unit  mass;  upper  case 
refers  to  total  quantities  in  a  storm.)   So  about  one  third  of  the  heat  goes 
into  expansion.   Using  the  1.3  x  1016  J  obtained  above  for  dH  and  p  =  5  x  104 

N  m~2,  we  get 

dV  =  7.5  x  1010  m3.  (3.5) 

From  the  equation  of  state  (with  pressure  constant)  and  (3.2)  we  obtain 

dV    dT    dh 


V     T     C  T 
P 


(3.6) 


A  typical  mixing  ratio  inside  a  severe  storm  is  up  to  10  gm  H2O  per  kilogram 
of  air,  giving  a  latent  heat  release  dh  of  about  2.5  x  104  J  kg~l  of  air, 
assuming  all  the  vapor  condenses.   Using  Cp  =  1005  J  kg-1  K"1,  we  get  for 
the  fractional  volume  and  temperature  changes, 

dV     dT    _  __  ,_  _N 

—  =  —  =  0.09,  (3.7) 

or  almost  10  percent.   Such  a  volume  change  occurs  over  the  storm  growth 
phase,  nominally  1000  s. 

The  spectral  content  of  the  acoustic  radiation  depends  on  the  exact  form 
of  V(t),  the  temporal  behavior  of  the  expanding  volume;  if  it  is  smooth,  the 
emissions  will  contain  mainly  the  frequency  characteristic  of  the  storm  growth 
time.   Higher  frequencies  would  be  generated  by  fluctuations  in  volume  change 
(caused,  for  example,  by  fluctuations  in  heat  release)  superimposed  on  the 
overall  storm  growth. 

I  now  estimate  an  upper  bound  for  the  acoustic  energy  contained  in  mono- 
pole  radiation  at  various  high  frequencies.   Let  Vm  be  the  total  volume  change 
(caused  by  latent  heat  release)  that  occurs  during  the  storm-growth  time  T. 
Suppose  volume  fluctuations  of  shorter  period  T1  are  superimposed  on  this 
growth,  but  impose  the  additional  condition  that  V  must  increase  monotonically 
with  time  (this  means  that  heat  is  not  periodically  fed  back  to  evaporate 
water  droplets;  we  will  discuss  this  possibility  separately  in  Sec.  7).   The 
amplitude  of  volume  fluctuations  at  wave  period  T'  cannot  be  much  larger  that 

T' 

=—  V  .   The  time  rate  of  change  of  volume  at  period  T'  is  thus  of  the  order 

T     j  m  T, 

of  7zt    sr~     V    ,    or  just   the   same   as   the  average  rate  of  change  of  volume.      With 
T*      T         m  J  to  ° 
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the  value  of  Vm  given  in  Eq.  (3.5)  and  a  1000- s  lifetime,  we  get  an  average 
volume  rate  of  change 

S  =  — -  a  7.5  x  107  m3  s"1  .  (3.8) 

Morse  and  Ingard  (1968,  Sec.  7.1)  give  the  following  formula  for  the  rms 
acoustic  radiated  power  II  from  a  monopole  source  of  strength  S  at  frequency  co 

n  =  g^_  |S  |2  ,  (3.9) 

4ttC  '  co '   '  J 

where  C  is  the  speed  of  sound  and  p  is  the  air  density.   Leaving  the  wave 
frequency  oj  variable  for  the  moment,  and  using  values  appropriate  to  the 
troposphere,  we  get  (in  mks  units) 

II  *  1.6  x  1012  w2  .  (3.10) 


If  all  the  energy  were  radiated  at  a  wave  period  of  30  s,  for  example, 
(oo  =  0.209  s_1)  II  =  7  x  1010  W.   This  is  at  least  two  orders  of  magnitude 
larger  than  the  infrasonic  wave  observations  require. 

At  the  present  time,  there  is  no  theoretical  or  observational  reason  to 
believe  that  such  coherent  higher  frequency  fluctuations  in  latent  heat  re- 
lease actually  occur.   All  we  have  shown  is  that  such  fluctuations  could 
easily  account  for  the  observed  wave  intensities.   Arguing  against  coherent 
fluctuations  in  latent  heat  release  over  such  a  large  volume  is  the  fact  that 
the  typical  storm  size  is  about  an  acoustic  wavelength  at  tens  of  seconds  wave 
period.   A  fluctuation  in  one  part  of  a  storm  would  change  phase  substantially 
before  its  effects  could  be  communicated  acoustically  across  the  storm. 


3.2  Dipole  Source 

Not  all  the  motions  inside  a  storm  are  caused  directly  by  the  release  of 
latent  heat  of  vaporization;  we  can  conceptually  separate  the  fundamentally 
expansive  motions,  which  we  represented  by  an  acoustic  monopole,  from  the 
cellular  or  (ideally)  toroidal  convection  pattern  caused  partly  by  the  static 
instability  of  the  atmosphere  and  reinforced  by  latent  heat  release.   I  now 
estimate  the  equivalent  acoustic  dipole  strength  using  measurements  of  thun- 
derstorm updraft  speed  and  size,  but  again,  the  spectral  content  is  unknown  on 
time  scales  shorter  than  the  storm's  lifetime. 
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There  are  several  practically  equivalent  ways  to  model  the  idealized 
convection  pattern  of  a  thunderstorm  that  would  permit  us  to  estimate  its 
dipole  radiation.   First,  we  could  assume  a  sphere,  or  bubble,  of  fluid  accel- 
erating vertically;  second,  we  could  assume  a  thin  horizontal  disk  with  vary- 
ing fluid  dipole  sources  distributed  over  its  surface;  or  third,  we  could  use 
a  horizontal  vortex  ring  whose  strength  varies  with  time.   All  three  methods 
yield  essentially  the  same  hydrodynamic  flow  field  and  the  same  acoustic 
radiation  for  a  given  spectral  content  of  the  source  variable.  I  use  the  first 
method  for  a  sample  calculation. 

Morse  and  Ingard  (1968,  Sec.  7.1)  give  a  formula  for  the  acoustic  power  II 
radiated  by  a  sphere  of  radius  a_  oscillating  up  and  down  with  a  frequency  CO 
and  maximum  speed  U: 

4 
n  =  pC(4Tra2)  ^—    U2  ,  (3.11) 

where  k  =  co/C.   Evaluating  the  constants  in  mks  units  for  the  troposphere 

n  =  3.2  x  10"8  a6oo4U2  Watts.  (3.12) 

-1  12 

For  a  =  2.5  km,  U  =  10  m  s   and  a  wave  period  of  30  s,  II  =  1. 5  x  10   W.  If 

we  take  the  estimate  of  the  average  broadband  acoustic  radiated  power  required 

by  the  infrasonic  measurements,  about  108  to  10$  W,  and  the  same  values  for  a 

and  co,  we  require  velocity  fluctuations  of  only  0.08  to  0.25  m  s-*  at  30  s 

period,  not  unreasonable  for  fluctuations  in  thunderstorm  updrafts.  Note, 

however,  that  the  model  requires  the  fluctuations  to  be  coherent  over  a  sphere 

of  radius  a.      If  we  decrease  the  value  of  a  by  a  factor  n,  the  required  U 

increases  by  a  factor  of  n^. 

The  remarks  at  the  end  of  the  last  section  can  be  paraphrased  here:  a 
dipole  source  appears  to  be  quantitatively  realistic  only  if  coherent  motions 
of  the  required  amplitude  exist  at  the  proper  frequencies.  Again,  supporting 
experimental  data  are  lacking. 


3. 3  Quadrupole  and  Higher-Order  Sources 

I  could  extend  the  discussions  of  the  preceding  sections  to  calculate  the 
radiation  from  higher-order  simple  sources;  however,  this  would  probably  serve 
no  useful  purpose,  because  no  identifiable  storm  processes  suggest  description 
by  such  sources,  and  because  the  same  problems  arise  with  respect  to  the 
energy  in  higher- frequency  fluctuations. 

Models  of  sound  radiation  by  turbulence  employ  a  random- quadrupole  repre- 
sentation of  fluctuations  in  Reynolds  stresses  in  fluids,  and  this  stochastic 
approach  will  be  examined  next. 
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4.   TURBULENT  SOURCES 

From  the  point  of  view  of  acoustic  radiation,  it  is  useful  to  divide 
fluid  motions  into  three  kinds:   (a)  steady,  or  time- invariant,  flows,  (b) 
nonsteady  flows  of  low  Reynolds  number,  whose  principal  dynamic  features  can 
be  described  "exactly,"  and  (c)  nonsteady  flows  of  high  Reynolds  number  that 
are  too  chaotic  to  be  described  deterministically.   That  steady  flow  radiates 
no  sound  is  obvious,  for  example,  from  energy  considerations,  and  it  is  there- 
fore of  no  direct  interest  to  us.   The  sections  on  simple  sources  and  on 
vortex  sound  contain  examples  of  nonsteady  flows  that  are  modeled  in  an  exact 
way,  and  whose  radiated  sound  is  exactly  determined  from  the  flow  field.   The 
third  kind  of  flow,  which  we  can  at  present  describe  only  in  terms  of  its 
statistically  averaged  properties,  is  called  "turbulence."  Turbulent  flows 
also  radiate  sound,  but  because  the  flow  is  described  only  statistically,  the 
sound  field,  too,  is  described  only  in  statistical  terms. 

The  airflow  inside  severe  thunderstorms  indeed  defies  exact  description 
and  so  qualifies  for  the  label  "turbulent";  yet  certain  gross  storm  features 
seem  to  fit  deterministic  models  for  convective  flow,  and  so  in  some  sense  can 
be  considered  "known."   (It  is,  in  fact,  a  steady  migration  of  storm  charac- 
teristics from  the  former  category  into  the  latter  that  measures  what  we  call 
an  increased  "understanding"  of  storms.)   Yet  it  appears  that  storm  airflow 
will  never  be  describable  in  detail,  and  that  the  ultimate  storm  model  will  be 
some  deterministic  description  of  storm-scale  flow,  upon  which  are  superim- 
posed smaller-scale  "turbulent"  fluctuations,  whose  influence  is  parameterized 
only  in  statistical  terms.   The  question  we  face  here  is  whether  the  deter- 
ministic part  or  the  random  part  of  storm  motions,  or  some  combination  of 
both,  is  required  to  explain  the  observed  acoustic  radiation.   To  answer  this 
question,  I  next  examine  the  predictions  of  a  turbulent- source  model.   A 
turbulent  description  of  the  severe-storm  radiation  mechanism  is  not  very 
rewarding  from  the  point  of  view  of  insight  into  storm  processes,  but  it  can 
be  used  to  estimate  average  radiated  power  and  the  average  acoustic  spectrum 
when  no  exact  description  is  possible. 

To  construct  a  simple  turbulent  description  of  severe-storm  air  motions, 
suppose  that  energy  is  supplied  to  the  storm's  motions  (by  convective  insta- 
bility) on  an  "outer  scale",  L,  equal  to  the  scale  of  the  convection  cell 
itself.   According  to  the  Kolmogorov  hypothesis  (Lumley  and  Panofsky,  1964), 
motions  on  smaller  scales  merely  partake  of  an  energy-cascade  process  that 
transfers  energy  to  smaller  and  smaller  scale  motions,  until  an  "inner  scale" 
is  reached  where  energy  can  be  efficiently  dissipated  by  molecular  viscosity. 
In  the  in-between,  or  so-called  inert ial  range  of  scales,  where  energy  is 
neither  supplied  to  nor  removed  from  the  fluid  motions  but  is  merely  trans- 
ferred by  means  of  nonlinear  interactions  between  scale  sizes,  the  statistical 
and  spectral  properties  of  the  turbulence  (and  its  radiation)  do  not  depend  on 
the  geometry  of  the  driving  mechanism. 

Turbulence  in  the  inertial  range  is  said  to  be  homogeneous  and  isotropic, 
properties  that  simplify  its  analysis  and  the  calculation  of  its  acoustic 
radiation.   From  dimensional  arguments,  it  is  possible  to  establish  that  the 
frequency  spectrum  of  the  (squared)  pressure  fluctuations  in  such  a  fluid 
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exhibits  an  f     behavior  (Lumley  and  Panofsky,  1964),  but  that  the  power 

spectral  density  of  its  acoustic  radiation  falls  off  with  frequency  as  f~7/2 
(Morse  and  Ingard,  1968,  Sec.  11.4;  Meecham  and  Ford,  1958).   If  L  is  the 
largest  scale  of  fluctuations  in  the  storm  we  expect  the  radiated  power  spec- 
trum to  level  off  at  a  wave  frequency  of  about  Ul/2ttL,  where  Ul  is  a  repre- 
sentative velocity  fluctuation  on  the  scale  L  and  then  to  fall  off  again  at 
lower  frequencies.  For  the  largest  thunderstorms,  Ul  -  30  m  s"1  and  L  -  10  km, 
which  give  a  wave  period  of  about  30  min.   So,  it  seems  reasonable  to  hypothe- 
size a  wave  spectrum  for  thunderstorm  radiation  that  peaks  in  the  vicinity  of 
20  to  30  min  (actually,  these  would  be  gravity  waves)  and  falls  off  at  lower 
and  higher  frequencies,  with  the  high-frequency  falloff  approaching  an  f-7/2 
rate.   Meecham  and  Ford  predict  a  low- frequency  spectral  density  proportional 
to  f4. 

For  frequencies  near  the  peak,  however,  the  assumptions  of  a  turbulent 
description  (i.e.,  that  the  flow  is  of  sufficient  spatial  and  temporal  extent 
to  permit  valid  averages  to  be  taken)  break  down.   There,  we  would  expect 
deterministic  descriptions  of  the  radiation  to  be  more  appropriate,  for  ex- 
ample, the  acoustic  monopole  and  dipole  models  described  in  Sec.  3,  applied  to 
scales  of  size  L.   A  complete  storm  radiation  model  might  then  consist  of  a 
turbulent  description  of  the  high-frequency  radiation,  matched  with  a  deter- 
ministic description  of  lower  frequency  radiation  based  on  known  storm  motions. 

To  estimate  the  amount  of  acoustic  power  that  could  be  radiated  by  tur- 
bulent motions  of  the  intensity  thought  to  occur  within  thunderstorms,  I  use 
an  approximate  formula  derived  by  Proudman  (1952),  based  on  Lighthill's  (1952) 
theory  of  acoustic  radiation  from  homogeneous,  isotropic  turbulence.  Light- 
hill's  hypothesis  was  that  the  sound  radiated  from  a  region  of  turbulent  fluid 
isolated  from  external  forces  is  equivalent  to  the  sound  radiated  by  a  random 
distribution  of  incoherent  acoustic  quadrupoles  whose  strength  is  related  to 
the  fluctuating  stress  field  in  the  fluid.   In  air,  on  the  scales  of  motion  of 
interest  to  us,  the  principal  form  of  fluid  stress  is  not  the  "real"  or  vis- 
cous stress  exerted  across  a  fluid  element  by  molecular  motions  normal  to  its 
boundaries;  instead,  distortions  of  fluid  elements  can  be  thought  of  as  caused 
by  Reynolds  stresses  -  that  is,  by  momentum  pu^  carried  across  the  boundaries 
of  a  fluid  element  by  fluid  motion  uj .   Because  such  stresses  in  a  fluid  tend 
to  distort  the  shapes  of  fluid  elements  and  cause  them  to  rotate,  one  can 
think  in  physical  terms  of  a  distribution  of  point  quadrupoles,  located  one 
per  average  eddy  volume.   The  acoustic  quadrupole  strength  is  proportional  to 
the  strength  of  the  fluctuating  part  of  the  Reynolds  stresses. 

The  final  result  of  Proudman' s  analysis  (see  also  Morse  and  Ingard,  1968, 
Sec.  11.4)  gives  the  total  acoustic  power  IT  radiated  by  a  volume  V  of  homo- 
geneous, isotropic  turbulence,  in  which  there  are  rms  Mach-number  fluctuations 
M  and  a  turbulent  dissipation  rate  e  per  unit  mass: 

n  =  38M5epV   .  (4.1) 
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The  value  of  £  indicates  the  intensity  of  turbulence  and  has  been  estimated 
inside  thunderstorms  from  Doppler  radar  measurements  (Strauch  et  al . ,  1975) 
to  be  between  0.1  and  1.0  m^s-^;  the  latter  value  applies  to  the  most  severe 
storms,  and  we  will  use  it.   (Independent  order- of- magnitude  estimates  of  e 
can  be  obtained  using  u$/l,   where  u  represents  velocity  fluctuations  on  any 
scale  I.)      We  assume  the  largest  reasonable  value  for  the  storm  volume,  V,  a 
10  km  cube,  even  though  we  recognize  that  the  assumed  values  of  £  and  M  are 
maximum  values  and  are  not  maintained  over  the  whole  storm  volume.  Using  a 
value  of  0.03  for  M  (about  10  m  s-*  velocity  fluctuations)  gives  II  =  9.2 
x  10*>  W;  a  value  of  0.1  for  M  (about  34  m  s""*  velocity  fluctuations)  gives 
n  =  3.8  x  108  W. 

As  mentioned  earlier,  this  power  should  be  distributed  (on  the  average) 
over  a  spectral  region  beginning  at  about  30-min  wave  period,  with  a  power 
spectral  density  that  falls  off  toward  higher  frequencies,  ultimately  ap- 
proaching an  f-7/2  power- law  slope.   The  total  effective  width  of  the  emission 
spectrum  seems  to  depend  on  the  shape  of  the  spectrum  in  the  transition  region 
between  the  frequencies  corresponding  to  the  outer  scale  and  those  scales  that 
can  be  considered  fully  in  the  inertial  regime.  If  the  transition  takes  place 
in  a  decade  or  less,  the  total  emission  spectrum  would  appear  quite  narrow; 
with  an  f~7/2  falloff,  more  than  80%  of  the  power  appearing  above  any  given 
frequency  is  concentrated  in  the  octave  immediately  above  that  frequency.   The 
power  density  at  the  higher  frequencies  of  interest  to  us  would  then  certainly 
be  inadequate.   On  the  other  hand,  if  the  transition  takes  several  decades, 
the  spectrum  could  still  be  relatively  flat  at  wave  periods  of  tens  of  sec- 
onds.  But  then  the  results  for  radiation  from  the  inertial  range  would  not 
apply.   I  see  no  way  to  predict  the  spectral  extent  of  the  transition  region 
on  theoretical  grounds  alone;  new  storm  observations,  such  as  those  now  being 
made  with  dual-Doppler  radars,  may  help  answer  this  question. 

There  are  so  many  uncertainties  associated  with  the  turbulent -source 
model  that  it  is  not  now  possible  to  accept  or  reject  it  conclusively  as  a 
viable  explanation  for  the  observed  emissions.   Order-of-magnitude  uncertain- 
ties exist  in  formulas  like  Equation  (4.1)  as  well  as  in  those  describing  the 
spectral  properties  of  the  emissions.   Large  uncertainties  also  exist  in  the 
storm  parameters  I  assumed,  because  they  have  not  yet  been  measured  accurately 
or  extensively.   Only  by  using  values  for  M,  £,  and  storm  volume  that  are 
probably  all  larger  than  those  actually  present,  even  in  the  most  severe 
storms,  can  we  obtain  values  for  the  acoustic  radiated  power  that  come  close 
to  the  measured  estimates. 

I  conclude  that  a  purely  turbulent  description  of  the  emission  mechanism 
is  probably  inadequate.   I  think  that  its  main  shortcoming  is  that  it  assumes 
a  random  distribution  of  incoherent  sources ,  whereas  there  appear  to  be  storm 
processes,  vortex  motion  for  example,  that  could  produce  coherent  fluctuations 
over  scales  substantially  larger  than  u/f,  where  f  corresponds  to  the  radiated 
frequency  and  u  is  the  fluctuation  velocity  at  that  frequency.   Infrasound 
waveforms  that  often  exhibit  many  apparently  coherent  wave  cycles  also  provide 
evidence  against  a  purely  stochastic  viewpoint. 
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ELECTRICAL  SOURCES 


5.1  Lightning 

One  obvious  source  of  sound  in  severe  storms  is  thunder,  the  acoustic 
remnant  of  the  shock  wave  generated  by  electrical  discharges  (lightning).   It 
could  reasonably  be  argued  that,  whereas  the  audible  portion  of  the  thunder 
spectrum  is  attenuated  by  ordinary  absorption  in  distances  of  a  few  tens  of 
kilometers  (as  is  observed) ,  an  infrasonic  portion  (if  one  exists)  would 
travel  unheard  with  low  loss  to  hundreds  or  thousands  of  kilometers.   I  thus 
ask  first  whether  the  acoustic  spectrum  of  thunder  extends  to  infrasonic 
frequencies  with  sufficient  intensity  to  explain  the  observations. 

Holmes  et  al.  (1971)  have  measured  the  mean  total  acoustic  energy  of  30 
cloud-ground  flashes  as  6.3  x  106  J,  and  Remillard  (1960)  gives  4  x  107  J  as 
the  acoustic  energy  radiated  by  a  typical  flash.   The  most  energetic  flash 
measured  by  Holmes  et  al.  had  an  estimated  total  acoustic  energy  of  17  x  106  J. 
These  estimates  apply  to  the  smaller,  isolated  thunderstorms  such  as  occur  in 
the  New  Mexico  area.   Somewhat  larger  values  may  apply  to  Great  Plains  frontal 
storms,  but  no  measurements  seem  to  be  available. 

How  much  of  this  energy  lies  in  the  infrasonic  portion  of  the  spectrum? 
Few  direct  measurements  of  the  infrasonic  spectrum  of  nearby  thunder  have  been 
attempted,  but  those  available  suggest  a  sharp  falloff  in  power  density  below 
audible  frequencies.   For  example,  Few  et  al.  (1967)  used  pressure  sensors 
with  essentially  flat  response  down  to  0.1  Hz  and  reported  a  sharp  falloff  in 
power  density  below  10  Hz.   This  and  other  similar  spectra  indicate  that  the 
total  amount  of  energy  below  about  1  Hz  is  several  orders  of  magnitude  less 
than  the  total  acoustic  energy  radiated.   This  would  put  the  total  acoustic 
energy  per  flash  below  1  Hz  at  no  more  than  about  103  or  104  J,  hardly  enough 
to  account  for  the  10$  to  10^  W  (J  s-l)  of  average  broadband  source  power 
required  by  infrasonic  observations  (Part  I).  Furthermore,  it  is  difficult  to 
reconcile  the  impulsive  nature  of  thunder  with  the  observations  of  nearly 
continuous  emissions,  even  from  nearby  storms,  where  propagation  effects 
cannot  disperse  a  pulse  train. 

I  am  aware  of  no  attempt  to  measure  the  acoustic  spectrum  of  thunder  that 
has  utilized  instruments  suitable  for  the  portion  of  the  spectrum  (below  0.1 
Hz)  that  the  severe-weather  waves  occupy.   In  particular,  the  effects  of  wind 
noise  increase  rapidly  at  wave  periods  of  tens  of  seconds,  and  rather  elab- 
orate noise-suppression  techniques  must  be  used.  Therefore,  the  question  of 
the  infrasonic  energy  content  of  thunder  must  be  considered  open,  although  it 
seems  unlikely  that  significant  infrasound  radiates  from  the  lightning  channel 
itself.   Another  mechanism  has  been  proposed,  however,  for  low- frequency 
acoustic  pulse  radiation  during  thunderstorm  electrification. 
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5.2  Electrostatic  Sound 

Discussion  is  continuing  in  the  literature  about  whether  the  sudden 
release  of  electromechanical  stresses  in  clouds  contributes  to  the  infrasonic 
spectrum  of  thunder  (Colgate  and  McKee,  1969;  Holmes  et  al.,  1971a,  1971b; 
Bhartendu,  1971;  Dessler,  1973).   In  Dessler's  model,  electromechanical 
stresses  build  up  as  charged  water  droplets  rearrange  themselves  inside  a 
cloud  to  maintain  a  balance  between  the  electrostatic  forces  (arising  from 
charge  separation)  and  the  dynamic  pressure  each  droplet  experiences  by  its 
motion  through  the  air.   When  the  electrostatic  field  suddenly  changes  after  a 
lightning  discharge,  each  droplet  adjusts  its  motion  in  response  to  the  sud- 
denly uncompensated  pressure  field  in  its  vicinity.   This  motion  radiates 
sound.   In  simple  calculations,  Dessler  considers  disk-  and  cylinder- shaped 
volumes  of  like-charged  droplets,  initially  at  rest  but  which  begin  to  move 
apart  by  mutual  repulsion.   The  resulting  dynamic  pressure  field  contains  a 
region  of  reduced  pressure  at  the  center  of  the  disk  (or  cylinder) .  In  perhaps 
oversimplified  terms,  the  moving  droplets  "sweep  some  air  along  with  them, 
thus  creating  a  disk-shaped  region  of  reduced  pressure."  In  these  examples,  a 
sudden  neutralization  of  the  charge  distribution  (by  a  lightning  discharge) 
causes  a  momentary  implosion  that  radiates  an  acoustic-rarefaction  pulse. 

In  an  actual  thunderstorm,  the  charge  distribution  and  the  dynamics  of 
droplet  motion  are  probably  drastically  modified  by  the  dynamical  processes 
that  produce  and  maintain  charge  separation,  processes  whose  nature  remain 
obscure.   Therefore,  the  particular  geometries  treated  by  Dessler  must  be 
regarded  as  idealized  and  possibly  unrealistic.   Furthermore,  Dessler  did  not 
calculate  the  acoustic  energy  radiated  by  the  imploding  volume;  his  estimates 
of  signal  amplitude  (0.5  to  50  ubar)  are  actually  estimates  of  the  pressure 
deficit  at  the  center  of  the  volume  and  not  the  amplitude  of  the  radiated 
acoustic  wave. 


5.2.1  A  Model  for  Radiation  From  an  Imploding  Sphere 

A  more  satisfactory  estimate  of  radiated  acoustic  energy  uses  Lamb's 
(1945)  result  for  the  pressure  disturbance  radiated  by  a  sphere  of  radius  a 
with  an  initial  overpressure  Ap  inside,  which  is  released  at  t  =  0.  The  loss 
of  realism  suffered  from  the  use  of  a  spherical  volume  (rather  than  the  disk- 
shaped  or  cylindrical  models  Dessler  uses)  seems  minor  except  possibly  as  it 
affects  the  acoustic  radiation  pattern. 

Lamb  shows  that  the  density  disturbance  p1  has  the  form 

p,  =  Ap  o^co     f  (51) 

where  Ap  is  the  excess  density  inside  the  sphere  at  time  t  =  0,  C  is  the 
speed  of  sound,  and  r  is  the  distance  from  the  center  of  the  sphere. 
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R 

Assuming  the  initial  state  is  isothermal  and  choosing  a  time  t  =  ~—  ,  at 

which  the  spherical  disturbance  has  reached  a  radius  R  »a,  we  obtain  the 
classical  N-shaped  pressure  disturbance  (Figure  3) : 

A       R 
pt  =  I^E  (i  _  _°  )      r  _a<r<R  +a 
r     2        r         o      o 

p'  =  0  ,     Ro-a>r>Ro+a  ,  (5.2) 

where  y  is  the  usual  specific  heat  ratio.   Evaluating  (5.2)  at  r  =  R  ±  a 
gives  the  peak  value  of  the  pressure  disturbance, 

p,    4-l-=^L      R»a.  (5.3) 

r  max    2R±a      2   R       o  v   ' 

o  o 

We  expect  the  real  initial  pressure  distribution  to  be  continuous  (i.e., 
smooth)  across  the  sphere  boundary,  which  would  round  off  the  sharp  edges  of 
the  N-shaped  disturbance  somewhat.   In  the  case  of  interest  to  us,  Ap  is 
negative,  so  the  initial  pressure  disturbance  would  be  a  rarefaction. 

To  obtain  the  total  energy  radiated,  Er,  we  integrate  the  energy  density 
p'2/p  C^  over  the  spherical  wavefront  of  radius  RQ  from  R0-a  to  R0+a: 

2 

E  «    °  A)2       /  °    (l--)2  dr.  (5.4) 

r   p  C2    2    R  -a       r 
o  o 

3  2 
The  value  of  the  integral  is  8a  /R  ,  to  second  order  in  a/R  ,  so  that 

_  8-rra3  ryAp  .2 

r  "  p    L  C  J 
o 

which  is  independent  of  R   (for  R  »a) ,  as  it  should  be. 

We  can  also  estimate  the  acoustic  radiation  efficiency,  n,  of  an  explod- 
ing (or  imploding)  sphere  by  comparing  the  radiated  energy  with  the  work  done 
in  an  adiabatic  expansion  of  the  sphere.   For  adiabatic  processes,  pv"Y  = 
constant,  where  v  is  the  specific  volume  of  the  gas.  Logarithmic  differentia- 
tion shows  that 

pdv  =  -  —  dp.  (5.5) 
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r=R, 


Figure  3.      The  -pressure  disturbance  p'  as  a  function  of  distance     r 
radiated  by  a  sphere  of  radius  a  and  overpressure   Ap.,  which  is  re- 
leased instantaneously.      For   Ap  negative 3    the  rarefaction  phase 
precedes  the  condensation  phase. 


The  lefthand  side  is  the  work  done,  and  if  the  fractional  changes  in  p  and  v 
are  small,  the  work  done  is 


w=/pdv--^Ap 


(5.6) 


per  unit  mass.   For  the  entire  volume  V, 


W  =  -  ^  AP> 


(5.7) 


where  the  minus  sign  indicates  that,  for  Ap  positive,  work  is  done  by  the 
expanding  volume  on  the  surrounding  gas, 
radiation  efficiency, 


The  ratio  of  E  to  W  gives  the 


n  =  2]L   ^E=2y2  AR, 


(5.8) 


since  C  =  yd  /p  , 
ro  o 
overpressure. 


(T    Ko       ro 
The  radiation  efficiency  is  nearly  equal  to  the  fractional 


Using  values  appropriate  to  the  electrostatic  mechanism,  I  estimate  the 
radiated  acoustic  energy,  the  acoustic  efficiency,  and  the  amplitude  of  the 
pressure  wave.   According  to  Dessler,  the  maximum  value  of  Ap  that  can  be 
generated  by  electrostatic  forces  is  about  5  N  m-2  (50  ubar) .  The  radius  of 
the  volume  is  taken  to  be  5  km,  about  the  same  size  Dessler  uses.  Other 
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constants  are  given  the  following  values:    C  =  340  m  s   ,  y   =  1.4,   pQ 
=  1.2  kg  m~3.  Then  Er  =  1.1  x  1(P  J.   The  radiation  efficiency  is  1.3  x  10"^, 
and  the  peak  pressure  amplitude  at  a  distance  of  50  km  is  0.35  N  m"^  or 
5.5  ubar.   These  values  probably  represent  the  maximum  possible  field  obtain- 
able with  this  mechanism.   At  greater  distances,  the  field  strength  falls  off 
more  slowly  than  the  1/r  dependence  given  by  (5.2),  because  of  the  ducting 
effects  of  atmospheric  refraction. 


5.2.2  Energy  Comparison 

By  comparing  this  theoretical  estimate  of  the  maximum  energy  that  could 
be  radiated  from  a  single  electrostatic  source  of  reasonable  size  with  esti- 
mates of  average  radiated  power  based  on  severe-storm-wave  observations  (Part 
1),  we  can  infer  a  minimum  rate  of  such  emissions  that  would  be  required  to 
explain  the  observations.   The  average  broadband  power  radiated  by  the  most 
powerful  emitting  storms  was  estimated  in  Part  I  as  10°  to  10^  W.   If  the 
electrostatic  source  model  can  emit  no  more  than  10y  J,  such  emissions  every  1 
to  10  s  would  be  required  to  yield  the  observed  average  power  flux.   If  the 
emission  bandwidth  of  the  actual  waves  is  confined  mainly  to  the  observational 
"windows"  (discussed  in  Part  I) ,  then  a  smaller  average  power  flux  and  thus  a 
lower  emission  rate  would  be  required.   On  the  other  hand,  if  the  size  of  the 
emitting  volume  or  the  pressure  deficit  inside  it  is  smaller  than  the  nearly 
maximum  allowable  values  postulated  in  this  estimate,  then  the  radiated  energy 
would  be  considerably  smaller,  and  a  much  higher  emission  rate  would  be  re- 
quired.  On  an  energy  basis  then,  I  conclude  that  the  electrostatic  mechanism 
is  marginally  satisfactory,  but  there  are  large  uncertainties  in  the  model, 
and  to  some  extent  in  the  observations,  that  make  agreement  highly  tentative. 


5.2.3  Spectral  Considerations 

To  supplement  the  energy  comparisons  just  made,  I  .examine  possible  spec- 
tral discrepancies  between  model  and  observation.   As  in  the  radiation  from 
the  lightning  stroke,  the  sound  of  electromechanical  relaxation  processes  is 
essentially  impulsive,  depending  on  the  lightning  discharge  to  set  it  off. 
But,  in  contrast  with  the  sound  of  thunder,  the  size  of  each  disturbance 
radiated  by  electromechanical  relaxation  puts  its  spectral  peak  in  the  infra- 
sonic  portion  of  the  spectrum.   As  Lamb  showed,  the  length  of  the  pulse  radi- 
ated from  a  sphere  under  pressure  is  equal  to  the  sphere's  diameter.   The 
dominant  wave  period  contained  in  the  N-shaped  impulse  is  thus  2  a/C.   In  the 
case  of  a  sphere  of  5-km  radius,  the  acoustic  pulse  length  would  be  about  30  s, 
precisely  in  the  spectral  region  of  the  observed  waves.   If  the  observed 
infrasound  is  composed  of  the  superposition  of  such  N-waves,  somewhat  smoothed 
and  dispersed  by  various  propagation  effects  (e.g.,  selective  absorption, 
ducting) ,  then  we  would  expect  pressure  fluctuations  with  rapidly  changing 
amplitude  and  spectral  content,  arriving  from  quickly  changing  directions.   A 
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difficulty  arises,  however,  with  emission  rates  greater  than  about  once  every 
ten  seconds,  because  the  spectral  content  of  the  emissions  would  then  reflect 
the  repetition  rate  instead  of  the  wavelength  determined  by  the  source  size. 


5.2.4  Objections  to  the  Mechanism 

The  most  serious  objection  to  be  raised  against  the  electrostatic  mech- 
anism concerns  the  pressure  waveform  observed  from  nearby  storms  (see  Figure 
4).  At  such  distances,  propagation  effects  have  negligible  influence  on  the 
emitted  waveform,  which  should  be  received  virtually  intact  and  should  reveal 
the  impulsiveness  and  intermittency  that  the  electrical  mechanism  predicts.  In 
fact,  infrasound  from  nearby  storms  reveals  no  such  effects  (see  also  Bowman 
and  Bedard,  1971). 


Boulder    Microborographs  June  5,  1973 
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Figure   4a.      Four  superimposed  microbarograph  traces  recorded  at  Boulder 3 
Colo. j   during  a  severe -weather  outbreak    (including  tornadoes)   in  Okla- 
homa during  1973.      The  four  traces  have  been  shifted  in  time  for  maxi- 
mum cross-correlation;    the  shifts  correspond  to  delays  appropriate   to 
a  plane  acoustic  wave  passing  over  the  sensor  array  at  340  m/s  from 
128  deg  azimuth.      This  acoustic  event  was  monitored  at  Boulder  for 
over  5  hours  and  was  detected  at  six  other  observatories  as  far  away 
as  College j   Alaska. 
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Another  objection  can  be  raised  on  the  observational  grounds  that  infra- 
sound and  electrical  activity  are  not  observed  to  be  closely  related.  If  the 
infrasound  radiation  mechanism  were  fundamentally  electrical,  we  would  expect 
at  least  two  kinds  of  observational  verification:   First,  there  should  be  some 
statistical  relation  between  the  geographic  and  temporal  distributions  of 
infrasound  sources  and  lightning  activity.   No  such  associations  are  evident. 
For  example,  a  statistical  maximum  in  lightning  activity  exists  in  Florida, 
yet  no  significant  long-distance  infrasound  has  been  associated  with  Florida 
storms.   (Some  ionospheric  infrasound  has  recently  been  observed  in  Florida  by 
Prasad  et  al.  (1975),  however.) 


Boulder    Microbarograpris     1,2,4     June  4, 1976 

Shifted    for  o  wave   from  71°  at  360  m/s 


Mountain  Daylight  Time 


Figure  4b.      Three  superimposed  miorobaragraph  traces  recorded  at  Boulder 3 
Colo. j    during  the  passage  of  a  tornado  about  27  km  NE  of  the  observa- 
tory.     The   traces  have  been  shifted  in  time  corresponding  to  a  plane 
acoustic  wave  from  71  deg  azimuth  at  360  m/s3   which  coincides  with  the 
direction  of  the   tornadic  storm.      Compare  the  characteristics  of  the 
waveform  and  its  spectral  content  with   the   longer  range  event  of  Fig. 4a. 
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Second,  observable  radio- frequency  electromagnetic  radiation  should 
accompany  individual  infrasound  radiation  events.   An  experiment  to  observe 
such  a  connection  has  recently  been  carried  out,  and  its  results  are  reported 
in  another  paper  (Beasley  et  al.,  1976),  but  it  will  suffice  to  say  here  that 
no  evidence  of  a  connection  has  been  found. 

In  summary,  some  features  of  the  electrostatic  emission  mechanism  seem  to 
satisfy  observational  requirements;  yet  others  remain  unexplained.  Because 
explanations  can  be  conceived  for  the  failures,  I  resist  discarding  this 
mechanism  completely,  and  label  it  only  "unlikely." 


5.3  Nonlinear  Pulse-Stretching 

I  examine  one  more  process  by  which  lightning  discharges  could  be  respon- 
sible for  severe  storm  infrasound  --  the  possibility  that  nonlinear  pulse 
stretching  in  the  upper  atmosphere  shifts  the  energy  in  the  lightning  shock  to 
the  infrasonic  portion  of  the  frequency  spectrum,  where  it  can  propagate  to 
great  distances  relatively  unattenuated. 

This  would  work  as  follows:   Some  of  the  acoustic  energy  radiated  as 
thunder  reaches  the  upper  atmosphere  (100  to  120  km)  before  returning  to  earth 
(Figure  5).   In  propagating  to  such  altitudes,  the  weak  shock  experiences 
ez/2H  amplitude  growth  until  propagation  becomes  nonlinear;  the  positive  phase 
of  the  shock  waveform  travels  faster  than  the  lagging  negative  phase,  and  the 
pulse  becomes  stretched.   When  the  pulse  is  refracted  back  to  earth,  its 
amplitude  is  diminished  again,  but  its  length  remains  stretched  so  that  its 
spectral  content  is  shifted  to  lower  frequencies  (see  Figure  6) . 


Figure  5.     Acoustic  raypaths  in  a  rea- 
listic temperature  profile  with  no 
winds.    The  source  is  at  10  km  height, 
and  rays  are  spaced  5  deg  in  eleva- 
tion angle.    The  shaded  region  shows 
rays  that  reflect  in  the  100  to 
150  km  height  range 3   and  that  could 
experience  nonlinear  pulse- stretch- 
ing. 
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5.3.1  Pulse-Lengthening  Estimate 

Daniels  (1969)  has  adapted  a  formula  of  Reed's  (1959)  to  calculate  the 
amount  of  pulse  stretching  suffered  by  upward- traveling  shock  waves  from 
explosions  and  from  supersonic  aircraft  (sonic  booms) .   Both  our  everyday 
experience  and  some  measurements  reveal  that  the  fractional  overpressures  (at 
comparable  distances)  are  about  the  same  for  thunder  and  for  sonic  booms 
(Daniels,  1969;  Jones  et  al . ,  1968).   Daniels  gives  a  fractional  overpressure 
of  2.36  x  10~4  at  a  reference  distance  of  17.1  km  for  a  small  supersonic- 
fighter  aircraft  and  5.6  x  10"^  for  a  larger  aircraft  (SST) .   Extrapolating 
the  curves  of  Jones  et  al .  (their  Figure  1)  to  the  same  reference  range  shows 
that  the  range  of  values  given  for  sonic  boom  overpressure  corresponds  to 
"typical"  lightning  flashes  of  2  to  30  x  10^  J  km~l .  Their  spectral  contents 
(derived  from  the  length  of  the  positive  pressure  phase  of  the  pulse)  are  also 
comparable.   Therefore,  within  the  accuracy  of  our  knowledge  of  both  phenom- 
ena, the  pulse-stretching  calculations  Daniels  did  should  apply  to  lightning 
shocks  as  well.   He  calculated  pulse- lengthening  factors  of  13  to  17  times  for 
waves  reflecting  at  110  km;  the  peaks  of  the  frequency  spectra  were  shifted  to 
0.4  and  0.13  Hz,  well  into  the  frequency  range  of  interest  to  us.   Multiple 
reflections  between  the  ground  and  thermosphere  could  cause  further  lengthen- 
ing (Figure  6) . 


5.3.2  Absorption  Estimate 

To  find  out  if  waves  traveling  such  a  path  could  actually  be  observed,  we 
must  also  estimate  the  acoustic  absorption  suffered. 

Lacking  adequate  theoretical  or  experimental  information  about  nonclassi- 
cal  absorption  mechanisms  (e.g.,  by  turbulence  or  molecular  relaxation)  at 
infrasonic  frequencies,  we  estimate  absorption  caused  only  by  viscous  damping. 
The  results  should  serve  as  a  lower  bound  on  the  total  absorption  from  all 
mechanisms.   Several  sources  (cf.  Reed,  1972;  Procunier  and  Sharp,  1971)  give 
the  classical  sound  amplitude  absorption  coefficient  a  in  the  form 

a  =  kf2/p  ,  (5.9) 

where  f  is  the  wave  frequency,  and  p  is  the  ambient  atmospheric  pressure;  k  is 
a  constant  whose  mean  value  in  mks  units  is  about  2  x  10"6.   The  plane-wave 
acoustic  amplitude  is  then  calculated  from 

A  =  A  exp  (Z/2H)  exp  -  /  a   ds ,  (5.10) 

0  S 

over  the  path  S,  where  H  is  the  pressure  scale  height.   The  path  S  begins  near 
the  ground  and  terminates  at  a  height  Z.   Formula  (5.10)  accounts  for  the 
competing  effects  of  exponential  amplitude  growth  and  viscous  attenuation,  but 
ignores  spherical  divergence,  which  is  of  secondary  importance. 
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Acoustic  ray-tracing  calculations  (Figure  3)  show  that  if  acoustic  rays 
are  launched  from  10  km  altitude,  rays  having  elevation  angles  3  (from  hori- 
zontal) between  about  30  and  60  deg  travel  to  heights  between  about  100  and 
120  km  before  returning  to  the  ground.   These  are  the  rays  that  could  con- 


ceivably experience  pulse- lengthening  effects, 
following  estimates. 


We  assume  3  =  45  in  the 


For  heights  to  130  km  or  so,  we  can  assume  that  pressure  p  follows  an 
exponential  law  of  the  form 


p(z)  =  p(o)  e 


-z/H 


(5.11) 


where  z  is  the  height  and  H  is  the  pressure  scale  height,  given  a  constant 
value  6.65  km  here;  p(o)  is  10->  N  m~2.   We  simplify  our  absorption  estimate 
by  assuming  that  the  rays  travel  in  straight- line  triangular  paths,  reflecting 
at  a  height  Z.   The  elevation  angle  3  is  then 


3  =  tan   (2  n  Z/D) 


(5.12) 


where  D  is  the  total  ground  range  and  n  is  the  number  of  ray  "hops"  between 
the  ground  and  the  thermosphere  base.   Substituting  (5.9)  and  (5.11)  into 
(5.10),  we  get  the  up- leg  attenuation 


Zn  ~  =  Z/2H  -  ^-,  /  ez/H  ds. 
A  p(o)  J 

o  r    s 


(5.13) 


Substituting  dz  =  ds  sin3  and  integrating  from  z  =  0  to  Z ,  we  get 


-   A    7/ou    kf2H      f   Z/H  .. 

Zn  -r—  =  Z/2H  -  —7-^ 7—0      (e   -1) 

A  p(o)  snip 


(5.14) 


The  competing  effects  of  exponential  growth  and  viscous  damping  cause  the 
wave  to  reach  a  maximum  amplitude  at  the  height  where  the  growth  rate  equals 
the  damping  rate.   This  height  increases  with  decreasing  wave  frequency.  By 


differentiating  (5. 14)  with  respect  to  z  we  find  Z 
amplitude: 


max' 


the  height  of  maximum 


z    =  H  Zn  P(o)  Si"3  . 


max 


2H  k  f 


(5.15) 


Table  5.1  gives  values  of  Z    for  some  frequencies  of  interest. 
&  max 
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Table  5.1     Z         Values 
max 


f(Hz)  Height  of 

Maximum  Amplitude 
(km) 

0.1  128 

1.0  98 

10.0  67 

100.0  37 


These  heights  also  happen  to  be  very  nearly  the  heights  where  damping  alone 
reduces  the  wave  energy  by  3  dB: 


Z3dB  =  ('96  "  13*3  ln   ^  km'  (5.16) 

where  f  is  in  Hz.   In  general,  the  integrated  ground- to-upper- atmosphere 
attenuation  (in  decibels)  caused  by  viscosity  alone  is 

a  =  1.63  x  10"6  f2  ez/H  dB,  (5.17) 


valid  up  to  heights  of  about  130  km. 

Taking  z  =  100  km, 

a  =  5.54  f2  dB,  (5.18) 

giving  about  11  dB  attenuation  per  hop  at  1  Hz,  only  0.11  dB  per  hop  at  0.1 
Hz,  but  over  1100  dB  per  hop  at  10  Hz.   (The  effects  of  exponential  amplitude 
growth  on  wave  amplitude  cancel  out  on  the  up  and  down  legs  of  each  hop.) 

Given  a  maximum  stretching  factor  of  about  20  and  the  excessive  absorp- 
tion of  higher  frequencies,  the  emitted  waves  must  be  at  1  Hz  or  lower  fre- 
quencies, if  thunder  is  to  explain  the  infrasonic  observations  at  10  to  40  s 
wave  periods.  Yet  measurements  of  thunder  spectra  reveal  very  little  energy 
below  1  Hz. 


5.3.3  Summary  of  Further  Objections  to  the  Mechanism 

A.   The  lightning  source  is  impulsive  in  nature  and  would  require  fre- 
quent and  repeated  strokes  for  several  hours  to  account  for  the  continuous 
waves  that  are  observed. 
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B.  Severe-weather  infrasound  is  seen  at  distances  from  the  apparent 
source  storms  that  are  too  short  (less  than  300  km)  to  permit  thermospheric 
reflection. 

C.  The  ionospheric  effects  have  periods  of  2  to  5  min  and  persist  for 
hours.   The  surface  and  ionospheric  effects  have  such  similar  phenomenologies 
that  a  common  source  mechanism  is  suspected.   Pulse  stretching  could  not 
explain  the  low  frequencies  and  continuous  nature  of  the  ionospheric  effects. 


6.   VORTEX  SOUND 

As  noted  earlier,  steady,  subsonic  flows  (including  steady  vortices)  emit 
no  sound.   But  in  the  atmosphere,  all  organized  vortex  motion  appears  to  be 
unsteady  in  some  sense.   All  vortices  are  unsteady  during  their  formation  and 
their  decay,  and  although  the  two  processes  are  conceptually  distinct,  they 
probably  overlap  during  the  lifetimes  of  most  atmospheric  vortices.  Further- 
more, the  flow  inside  what  may  appear  to  be  a  steady  vortex  often  reveals 
smaller  scale  undulations  or  vortices,  or  it  may  contain  fully  developed 
turbulence.   These  observations  are  consistent  with  many  theoretical  analyses 
(cf .  Lord  Kelvin,  1910) ,  which  demonstrate  that  vortex  flow  is  unstable  with 
respect  to  many  kinds  of  radial  and  axial  perturbations.  Therefore,  we  would 
expect  most  atmospheric  vortices  to  radiate  acoustic  energy;  furthermore,  we 
would  expect  their  rotation  and  cylindrical  geometry  to  favor  the  sustained 
generation  of  particular  wave  frequencies. 

There  is  ample  evidence,  based  mostly  on  observations  of  motions  on 
scales  that  produce  audible  sounds,  that  vortices  do  radiate.   The  aeolian 
tone  (cf.  Powell,  1964;  Richardson,  1924)  is  the  classical  example.   There, 
vortices  that  form  in  the  wake  of  a  cylindrical  obstacle  placed  in  an  airflow 
are  shed  at  a  characteristic  rate  and  generate  sound  at  the  shedding  fre- 
quency.  Some  musical  instruments  are  based  on  a  phenomenon  known  as  the  "edge 
tone"  (Powell,  1961),  in  which  sound  is  radiated  by  the  vortex  pairs  that  form 
periodically  when  an  airstream  impinges  on  a  wedge-shaped  obstacle.  Chanaud 
(1969)  studied  the  aerodynamic  sound  radiated  by  a  rotating  disk  in  air  and 
found  that  periodic  vortex  motion  around  the  disk's  periphery  was  responsible 
for  some  of  the  lowest- frequency  radiated  sound.   The  same  sort  of  peripheral 
turbulence  (or  perhaps  more  orderly  peripheral  vortices)  could  explain  reports 
of  "musical  sounds"  emitted  by  tornadoes  (Abdullah,  1966) .   Dynamically  simi- 
lar motions,  but  on  larger  scales,  could  radiate  infrasound  that  would  go 
mostly  unobserved.   I  now  examine  how  mesoscale  vortices  form  and  dissipate  in 
the  atmosphere,  with  specific  attention  to  those  associated  with  convective 
storms.   Then  I  estimate  how  much  radiation  could  be  expected  from  two  simple 
examples  of  unsteady  vortex  flow. 
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6.1  Vortex  Formation 

Apart  from  mesoscale  vortices  that  are  mechanically  produced,  for  ex- 
ample, in  the  flow  past  mountains  and  islands  (Chopra,  1973,  and  references 
therein),  atmospheric  vortices  all  involve  the  coupling  of  ambient  vorticity 
with  strong  buoyant  convection  (Morton,  1966) .   But  convection-associated 
horizontal  convergence,  operating  on  the  background  vorticity  due  to  the 
earth's  rotation  alone,  is  insufficient  to  explain  the  observed  intensity  and 
development  times  of  the  mesoscale  vortices,  especially  tornadoes,  associated 
with  thunderstorms.   Such  vortices  must  form  as  a  result  of  vigorous  convec- 
tion in  a  neighborhood  with  a  richer  distribution  of  vorticity  than  that  due 
to  the  earth's  rotation  alone  (about  10~4  s-1  at  midlatitudes) .   In  the  case 
of  squall- line  thunderstorms,  additional  background  vorticity  is  available 
from  the  synoptic  situation.   In  the  vicinity  of  a  front,  strong  cyclonic 
circulation  is  often  present  on  a  synoptic  scale,  and  can  contribute  about  an 
additional  10~4  s~l  vorticity.   Furthermore,  frontal  systems  are  likely  to  be 
rich  in  the  potential  vorticity  that  arises  whenever  the  constant-density 
surfaces  are  inclined  to  the  constant-pressure  surfaces  (producing  solenoids) . 

The  horizontal  vorticity  associated  with  vertical  wind  shear  in  the 

-2_i 
boundary  layer  is  generally  around  10   s  x.   Convection  can  draw  up  and 

stretch  loops  of  formerly  horizontal  vortex  tubes  (the  so-called  "tilting 

term"  in  the  vorticity-tendency  equation) ,  causing  an  intensified,  erect 

vortex  pair  with  opposite  rotation  senses  (Davies-Jones  and  Kessler,  1974) . 

Such  pairs  are  now  commonly  seen  inside  storms  with  Doppler  radars  (cf.  Krop- 

fli  and  Miller,  1976).   Similar  vortex-pair  formation  has  been  observed  by 

Idso  (1975)  in  the  airflow  over  hills. 

Another  kind  of  unsteady  vortex  flow  is  favored  during  strong  convection 
when  certain  relationships  between  the  horizontal  inflow  angle  and  the  "con- 
figuration ratio"  are  satisfied.   (The  configuration  ratio,  or  aspect  ratio, 
is  the  ratio  of  the  diameter  of  the  zone  of  convection  to  the  height  of  the 
inflow  layer.)   In  laboratory  simulations  of  tornado- like  vortices,  Ward 
(1972)  demonstrated  conditions  that  favor  the  production  of  multiple  vortices 
rather  than  a  single  one: 

"In  the  laboratory,  more  than  one  vortex  can  be  produced 
only  when  the  configuration  ratio  is  greater  than  unity. 
When  the  ratio  is  4,  and  the  angle  of  inflow  (measured  from 
a  radial)  0  is  very  small,  2  to  3  deg,  a  single  vortex  forms; 
as  9  is  increased,  the  diameter  of  the  core  increases  pro- 
portionately.  The  scale  of  the  turbulence  also  increases  with 
the  inflow  angle,  and  when  the  latter  reaches  about  30  deg, 
a  vortex  pair  develops.   They  are  located  on  opposite  sides 
of  the  parent  vortex  near  the  radius  of  maximum  tangential 
speed,  and  the  pair  rotates  around  the  central  axis  at  about 
half  that  speed.   Once  the  pair  is  formed,  the  flow  becomes 
somewhat  more  stable,  and  the  cores  are  well  defined.   This 
mode  will  continue  when  the  inflow  angle  is  decreased  several 
degrees  below  that  required  for  pair  formation. 
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After  a  pair  is  formed,  a  third  vortex  will  form,  and  then 
a  fourth,  as  the  inflow  angle  is  further  increased.   In 
general,  the  greater  the  number  of  vortices  in  the  system, 
the  smaller  their  individual  core  diameters. 

The  critical  inflow  angle,  0C,  at  which  multiple  vortices 
form  is  very  sensitive  to  the  geometry  of  the  convection. 
A  vortex  pair  forms  at  an  inflow  angle  of  about  30  deg 
when  the  configuration  ratio  is  4;  when  that  ratio  approches 
unity,  the  angle  required  for  pair  formation  increases  to 
about  75  deg.   When  the  depth  of  the  convergence  layer  exceeds 
the  diameter  of  the  convection,  no  multiple  vortices  have  been 
observed." 

Davies-Jones  (1973)  reinterprets  Ward's  results,  stating  that  the  experi- 
mental transition  from  a  single  wide  vortex  to  a  pair  of  concentrated  vortices 
occurs  at  a  critical  "swirl  ratio,"  rather  than  a  critical  inflow  angle.  The 
swirl  ratio  is  the  ratio  of  maximum  tangential  speed  to  the  mean  updraft 
speed.   Evidently,  a  system  of  multiple  corotating  vortices  often  represents  a 
stable  configuration  under  some  ambient  conditions  in  which  a  single  vortex 
would  be  unstable.   Next,  we  find  that  similar  multivortex  systems  can  form 
during  the  breakdown  of  a  single  steady  vortex. 


6.2  Vortex  Breakdown  and  Dissipation 

A  simple  calculation  shows  that  storm-scale  atmospheric  vortices  (in- 
cluding tornadoes)  would  have  lifetimes  of  more  than  a  year  if  molecular 
viscosity  were  the  sole  dissipating  agent.   To  get  decay  times  compatible  with 
those  observed  (tens  of  minutes) ,  we  need  viscosity  coefficients  about  10^ 
greater  than  reasonable  values  for  the  molecular  viscosity  of  air;  that  is, 
about  1  to  10  m^  s~l  is  required.   If  we  attribute  this  damping  process  to 
eddy  viscosity,  we  can  estimate  the  eddy  size  le   required  by  ve/ue,  where  ve 
is  the  eddy-viscosity  coefficient  and  ue  is  a  typical  eddy-associated  air 
velocity.   For  realistic  velocities,  we  get  an  eddy  size  of  a  few  meters.   We 
can  thus  picture  a  vortex  surrounded  by  a  mixing  region,  perhaps  a  few  meters 
thick,  through  which  momentum  is  transmitted  to  the  environment  by  eddies  of 
various  sizes. 

There  is  now  substantial  evidence  that,  in  many  cases,  the  early  stages 
of  vortex  breakdown  and  dissipation  are  characterized  not  by  a  broad  spectrum 
of  eddy  sizes,  or  "turbulence,"  but  rather  by  oscillatory  and  quasi-periodic 
motions  resulting  from  various  kinds  of  vortex  instability  (Lord  Kelvin, 
1910).   In  fact,  several  kinds  of  observations  suggest  that  certain  peripheral 
wave  numbers  are  preferred  and  that  vortex  breakdown  into  a  small  number  of 
well-defined  smaller  vortices  may  be  common  in  the  atmosphere. 

Fujita  (1970,  1971)  has  shown  evidence  for  the  existence  of  "suction 
vortices"  around  the  peripheries  of  tornadoes  and  dust  devils.   These  suction 
vortices  appear  as  a  ring  of  smaller,  intense  vortices  that  rotate  about  the 
axis  of  a  weaker  parent  vortex.   Weske  and  Rankin  (1963)  demonstrated  in  the 
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laboratory  the  existence  of  well-defined  periodic  vortices  spaced  around  the 
periphery"  of  a  rotating  disk  in  water.   (A  solid,  rotating  disk  can  be  thought 
of  as  representing  the  viscous  core  of  a  parent  vortex.)  These  vortices  seem 
to  be  the  cylindrical  analogs  of  the  boundary- layer  eddies  that  form  in  the 
flow  over  a  plane,  rough  surface,  but  with  the  additional  requirement  of  an 
integer  number  of  eddies  around  the  disk's  periphery. 

The  fact  that  many  tornadoes  leave  sequential  ground  tracks  resembling 
cycloidal  orbits  is  interpreted  as  strong  evidence  that  many  tornadoes  are 
actually  sub-vortices  orbiting  inside  larger-scale  "tornado  cyclones"  (Agee 
et  al . ,  1976).   Earlier,  Darkow  (1971)  had  presented  a  statistical  analysis  of 
the  periodicity  of  tornado  production  in  selected  tornado  "families",  which 
suggested  a  mean  interval  of  45  min  between  successive  tornado  touchdowns. 

A  detailed  case  study  of  a  multifunnel  tornado  by  Agee  et  al .  (1975) 
revealed  a  parent  tornado  vortex  with  a  300  m  diameter,  around  which  rotated 
as  many  as  seven  secondary  vortices,  some  of  which  became  tornadoes.  These 
peripheral  vortices  were  observed  to  grow  or  dissipate,  depending  on  which 
side  of  the  parent  vortex  they  were  on  (with  respect  to  the  direction  they 
travel  as  a  whole),  and  they  had  30-  to  60-m  diameters.  The  tangential  velo- 
city of  the  subsidiary  vortices  about  a  common  center  varied  with  position  in 
the  parent  vortex  and  ranged  between  25  and  55  m  s   ,  corresponding  to  ground 
speeds  of  up  to  75  m  s~l.   The  secondary  vortices  thus  tend  to  bunch  up  near 
the  trailing  edge  of  the  parent  vortex. 

We  estimated  a  frequency  associated  with  the  secondary-vortex  motion  by 
examining  Figure  14  of  Agee  et  al .  (1975)  and  estimating  the  times  between  the 
appearance  of  successive  vortices  at  the  same  location  along  the  periphery  of 
the  parent.   Three  time  intervals  between  the  passage  of  successive  vortices 
past  the  9  o'clock  position  on  the  parent  were  10.3,  8.7,  and  3.2  s.   If 
acoustic  waves  were  radiated  from  the  system,  the  spectrum  of  the  radiation 
would  be  in  the  same  period  range  (see  Sec.  6.3.1).  Unfortunately,  no  micro- 
barographs  were  suitably  placed  to  detect  any  such  radiation  from  this  storm. 

A  picture  is  now  emerging  (Agee  et  al . ,  1976)  of  a  complex  hierarchy  of 
vortices  and  sub- vortices  associated  with  some  tornadic  storms.  The  hierarchy 
extends  through:   (a)  the  "tornado  cyclone",  perhaps  5  km  in  diameter,  which 
probably  coincides  with  "hook  echoes",  as  well  as  cyclonic  features  observed 
with  Doppler  radars  and  may  be  exhibited  visually  in  the  "wall"  or  "pedestal" 
clouds  often  accompanying  tornadoes;   (b)   a  "mini-tornado-cyclone",  perhaps 
several  hundred  meters  in  diameter,  corresponding  to  the  "parent"  of  multiple 
tornadoes;   (c)   tornadoes  themselves,  up  to  a  few  tens  of  meters  in  diameter; 
to   (d)   the  sub-tornado- scale  "suction  vortices",  up  to  a  few  meters  in 
diameter,  proposed  by  Fujita,  which  are  probably  responsible  for  the  most 
intense  damage.   Figure  7  illustrates  such  a  vortex  hierarchy. 

Such  a  picture  suggests  the  importance  of  the  role  played  by  dynamic  in- 
stabilities, as  opposed  to  turbulence,  in  the  decay  of  thunderstorm  vortices. 
Theoretically,  instabilities  in  cylindrical  flows  can  take  many  forms,  from  a 
simple  two-dimensional  roll-up  to  complex  three-dimensional  (e.g.,  helical) 
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Figure   7.     A  schematic  of  a  possible  hierarchy  of  multiple- 
vortex  systems  for  the  various  scales  of  motion  associated 
with  a  tornado-producing  thunderstorm   (after  Agee  et  al.s 
1976). 


structures  (cf.  Davy  et  aT.  ,  1968).   For  example,  a  vortex  is  unstable  with 
respect  to  radial  perturbations  if,  within  some  annulus ,  angular  momentum 
decreases  with  increasing  radius  from  the  axis  of  rotation  (Rayleigh's  cri- 
tierion)  (Yih,  1969;  Chandrasekhar,  1961)  .   Other  instabilities  resemble  the 
Kelvin-Helmholtz  instabilities  in  plane-stratified  shear  flow,  with  centri- 
fugal force  playing  the  role  of  gravity.   The  number  and  sizes  of  the  smaller 
vortices  into  which  a  larger  vortex  breaks  down  are  governed  by  the  intensity 
and  thickness  of  a  high-shear  annulus  and  its  circumference. 

The  same  laboratory  simulations  of  tornado- like  vortices  cited  earlier 
(Ward,  1972)  have  verified  that  vortex  breakdown  can  be  accompanied  by  the 
formation  of  a  few  vortices  of  next-smaller  scale: 

"In  the  experimental  apparatus,  when  (the  convection  diameter) 
approaches  a  critically  large  value,  a  single  vortex  be- 
comes unstable,  with  accompanying  large  scale  turbulence 
and  large  fluctuations  in  core  diameter.   It  appears  that  this 
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instability  is  the  result  of  the  presence  of  inward  directed 
forces  associated  with  radial  (inward)  momentum  flux  sufficient 
to  produce  a  vortex  core  whose  central  pressure  deficit  require- 
ment is  greater  than  can  be  maintained  in  the  presence  of  sub- 
sidence in  the  large  core.   In  this  condition  of  instability, 
when  one  or  more  vortices  of  smaller  horizontal  scale  form 
in  the  annulus  of  [Figure  8],  they  tend  to  be  stable,  since 
their  smaller  cores  can  maintain  a  greater  pressure  deficit". 


For  our  purposes,  it  is  not  so  important  to  understand  the  details  of  the 
instabilities  that  may  cause  vortex  breakdown  as  it  is  to  recognize  that  the 
breakdown  process  can  often  be  described  in  terms  of  relatively  orderly  and 
long-lived  subsidiary  flows. 


Figure  8.     A  sketch  of  the  surface 
flow  pattern  observed  by  Ward  in 
laboratory  modeling  of  vortex-pair 
generation. 


6.3  Acoustic  Radiation  by  Unsteady  Vorticity 

Powell  (1964)  developed  a  "theory  of  vortex  sound,"  in  which  the  sound 
radiation  from  a  low-speed  compressible  flow  can  be  estimated  when  only  the 
spatial  and  temporal  distribution  of  vorticity  is  known.  This  was  an  impor- 
tant advance  from  our  point  of  view;  we  can  now  compute  the  acoustic  field  in 
terms  of  observables  that  do  not  involve  the  fluid's  compressibility  or  other 
difficult-to-estimate  quantities  like  the  distribution  of  Reynolds  stresses  in 
the  fluid.   Some  features  of  the  airflow  in  and  around  storms  (such  as  those 
just  discussed)  lend  themselves  to  description  in  terms  of  simple  vorticity 
distributions. 
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Powell's  general  expression  for  the  far- field  acoustic  velocity  pertur- 
bation u(x)  at  a  distant  observation  point  x  is 

u(x)  =  -  — £-=■    /  y  L"  dV(y)*  .  (6.1) 

4ttx(T   V   X  X 
o 

The  quantity  L^  is  the  x  component  of  the  second  time  derivative  of  £  x  u, 
where  £  is  the  vorticity  and  u  is  the  velocity  field  in  the  source  volume 
V0;  yx~is  the  x  component  of  the  spatial  source- integration  variably  y;  x  is 
a  unit  vector  in  the  x  direction,  and  the  asterisk  indicates  the  usual  time 
retardation,  the  only~effect  of  the  fluid's  compressibility  explicitly  consi- 
dered.  This  approach  becomes  intuitively  palatable  if  the  analogy  with 
electromagnetic  fields  is  considered:   if  Z>   represents  an  electric  current 
distribution,  then  u  represents  its  magnetic  field. 

Powell  applied  his  result  to  the  radiation  from  a  corotating  vortex  pair 
and  to  the  aeolian-tone  radiation  by  vortices  shed  from  obstacles.   I  apply 
his  results  to  storm- scale  airflow  and  estimate  the  properties  of  the  radi- 
ated infrasound. 


6.3.1  Radiation  From  a  Corotating  Vortex  Pair 

When  two  or  more  vortices  form  near  each  other  they  interact,  roughly 
speaking,  as  though  each  were  a  passive  tracer  in  the  velocity  field  of  the 
other.   A  counterrotating  vortex  pair,  for  example,  undergoes  simple  trans- 
lation in  a  direction  perpendicular  to  the  plane  containing  their  axes,  and  a 
corotating  vortex  pair  rotates  about  a  common  axis  in  the  same  direction  that 
each  individual  vortex  rotates.   The  translating  pair  is  of  no  direct  inter- 
est acoustically,  but  the  rotating  pair  does  radiate. 

The  hydrodynamic  velocity  field  of  the  corotating  pair  is  qualitatively 
illustrated  in  the  streamline  pattern  of  Figure  9.   It  is  apparent  that,  as 
this  noncircular  streamline  pattern  rotates,  the  external  fluid  experiences  a 
radial  force  whose  amplitude  passes  through  two  cycles  for  each  single  rota- 
tion of  the  vortex  pair.   Acoustically,  then,  the  rotating  streamlines  pro- 
duce the  same  effect  on  the  external  fluid  as  a  rotating  lateral  quadrupole, 
as  shown  in  Figure  9.   Powell  (1964)  has  calculated  the  acoustic  radiated 
power,  II,  starting  with  Eq  (6.1),  and  using  the  following  model. 

In  a  fluid  of  density  p  and  sound  speed  C,  let  two  parallel  line  vor- 
tices of  axial  length  I,  equal  strength  T  and  separated  a  distance  2r  spin 
about  an  axis  midway  between  them  (Figure  10).  The  radiated  power  is, 

n  =  i|_  p-  (JL)8  i2.  (6.2) 
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Figure  9.      The  streamline  pattern  calculated  for  a  pair  of  corotating  line 
vortices  of  equal  strength   (left)  }   and  the  resulting  instantaneous 3     rota- 
ting acoustic  radiation  pattern   (right) 3   which  is  equivalent  to  that  of  a 
rotating  acoustic  quadrupole. 

The  quantity  r/(4fTr)  is  the  velocity  u  at  which  the  vortices  spin  about  the 
common  axis,  and  therefore  their  angular  spin  rate  is  oj  =  r/(47rr^).  Because 
the  spinning  vortices  generate  a  hydrodynamic  (and  also  acoustic)  field  like 
a  fluid  quadrupole,  the  radiated  frequency  is  twice  the  rotation  frequency: 
f  =  (jo/it  =  r/47T^r2  =  u/frr. 


By  comparing  the  power  and  frequency  given  by  the  spinning- vortex  model 
with  those  associated  with  infrasound  observations,  we  can  estimate  the 
required  vortex  parameters.   First,  we  see  that,  independent  of  frequency, 


u  = 


•14 


n1/8 

Kl/8zl/4 


(6.3) 


where  K  =  jf—  and  *-=■  =  8.97  x  10    for  standard  conditions  (MKS  units) 


Figure  10.     Geometry  and  notation  for  calculating  the  acoustic  radiation 
from  a  corotating  vortex-pair.     Symbols  are  defined  in  the  text. 
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Table  6. 1     Required  Vortex  Corotation  Speeds 


1=1   km  Z=3  km 

II (Watts)  u  (m  s-1)  u  (m  s-1) 


107  57  43 

108  76  57 

109  101  76 


7 
Table  6.1  gives,  for  I  =   1  km  and  3  km,  values  of  u  required  to  radiate  10  , 

10^  and  10§  W.   Because  this  radiation  mechanism  is  fundamentally  monochro- 
matic, we  can  compare  these  results  with  the  3  x  10^  w  source  power  estimated 
by  Georges  (1973)  from  both  the  ionospheric  and  surface  pressure  measurements. 
(The  10°  to  10^  W  estimated  source  power  we  used,  for  example,  in  the  turbu- 
lent source  calculations,  is  an  estimate  of  the  total  broadband  radiated 
power. )  A  3  x  10'  W  radiated  power  would  require  eddy  spin  speeds  of  about 
50  m  s--*-,  a  little  high,  but  not  in  serious  disagreement  with  motions  observed 
in  the  most  severe  storms  (Brown  et  al.,  1971;  Agee  et  al.,  1975). 

Since  f  =  u/irr,  we  can  calculate  the  value  of  r  for  the  previously  tabu- 
lated values  of  u  and  for  two  wave  periods,  30  and  200  s  (representing  the  ob- 
served wave  periods  on  the  ground  and  in  the  ionosphere).  Table  6.2  shows 
values  of  r  corresponding  to  the  values  of  u  given  in  Table  6.1. 


Table  6.  2     Required  Vortex  Separation 


u  (m  s  ) 


43 

57 

76 

101 


T  =   30   s 

T 

=  200  s 

r   (km) 

r 

(km) 

0.41 

2.74 

0.54 

3.62 

0.73 

4.87 

0.96 

6.42 
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Thus,  the  surface  observations  require  an  eddy  separation  2r  of  about  a 
kilometer,  and  the  ionospheric  observations  about  6  km.   Eddies  of  the  larger 
scale  are  now  visible  inside  thunderstorms  observed  with  dual-Doppler  radars 
(Kropfli  and  Miller,  1975,  1976),  but  structural  details  of  the  eddies  and 
their  time  variations  are  still  beyond  the  resolution  of  current  remote- 
sensing  instruments.   Count errot at ing  vortex  pairs  have  been  seen  with  Doppler 
radars  inside  some  thunderstorms;  they  have  diameters  of  up  to  5  km  and 
extend  in  height  through  most  of  the  storm.   The  required  spin  rate  for  two 
vortices  equals  one  half  the  acoustic  radiated  frequency,  but  would  be  cor- 
respondingly smaller  for  more  than  two  vortices. 

Using  the  circulation  theorem,  we  can  get  a  rough  indication  of  the 
average  vorticity  corresponding  to  each  case: 

C   =  ~  •  (6.4) 

av   r  v   ' 

For  example,  for  the  case  u  =  43  m  s   and  r  =  2.74  km,  £   =  0.03  s   ; 

for  the  case  u  =  101  m  s   and  r  =  0.96  km,  £   =  0.21  s   .   For  comparison, 

the  vorticity  in  the  core  of  a  typical  tornado  has  been  estimated  at  about 
1  s~l  (Davies- Jones  and  Kessler,  1974)  ,  and  the  maximum  vorticity  measured 
in  a  Colorado  thunderstorm  with  Doppler  radars  is  about  0.01  s~l  (Kropfli 
and  Miller,  1976). 

The  radiated  power  depends  only  on  the  speed  that  the  vortices  spin 
around  each  other,  but  as  the  system  becomes  larger  while  keeping  the  radi- 
ated power  constant,  the  angular  spin  rate,  and  thus  the  radiated  frequency, 
is  lowered.   In  other  words,  the  wavelength  of  the  radiation  is  related  to 
the  vortex-pair  geometry  by 


—  =  -  .  (6.5) 

irr   u 

Equation  (6.1)  tells  us  that,  in  a  turbulent  flow,  where  eddies  of  all 
sizes  are  present,  the  most  powerful  acoustic  radiation  tends  to  be  from  the 
largest  eddies.   Because  of  the  separation  factor  yx,  it  is  the  moment  of 
vorticity  that  generates  sound.   In  the  case  of  thunderstorms,  then,  the 
motion  of  storm-scale  eddies  would  cause  the  most  intense  radiation.   This  is 
consistent  with  the  prediction  of  a  turbulent  description  of  the  storm's 
motions. 

The  rotating-pair  model  is  only  one  example  of  many  kinds  of  unsteady 
vortex  motion.   Because  vortex  motion  is  unstable  with  respect  to  almost  any 
kind  of  perturbation,  a  corresponding  variety  of  possible  acoustic-radiation 
mechanisms  exists.   Radial  oscillations  and  helical  instabilities  could  be 
modeled  exactly  and  their  radiation  predicted,  but  we  have  no  observations  to 
indicate  which  model  to  choose.   In  any  event,  we  would  expect  results  quan- 
titatively similar  to  those  obtained  with  the  corotating  pair.   Other 
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multiple-vortex  systems  can  assume  complicated  rotating  configurations, 
especially  if  many  vortices  of  different  strength  and  orientation  interact. 
Their  acoustic  radiation  would  be  correspondingly  complicated,  and  a  turbu- 
lent model  would  be  required  to  represent  all  but  the  simplest  configura- 
tions. 


6.3.2  Aeolian-Tone  Radiation 

Assume  that  some  part  of  a  convective  storm,  probably  the  intense  up- 
draft  core,  forms  an  obstacle  and  diverts  the  middle- level  wind  flow.  This 
assumption  is  supported  by  the  observations  by  Fujita  and  Arnold  (1963)  and 
by  Jessup  (1972) ,  who  found  blocking  of  the  flow  on  the  upwind  side  of  storms 
and  velocity  maxima  on  their  flanks. 

Next  assume  that  suitable  conditions  frequently  exist  for  the  generation 
and  shedding  of  coherent  vortices  from  the  storm-obstacle,  in  a  manner  dyna- 
mically similar  to  the  production  of  Karman  vortex  streets  on  a  laboratory 
scale.   This  assumption  can  be  justified  on  observational  grounds  from  space 
observations  of  trailing  vortices  in  the  lee  of  islands,  made  visible  by 
cloud  patterns  (cf.  Chopra,  1973;  Zimmerman,  1969,  and  the  references  they 
cite).   The  radar  observations  by  Jessup  (1972),  and  by  those  he  cites, 
appear  to  show  similar  vortices  in  the  lee  of  some  thunderstorms. 

The  apparent  dynamic  similitude  of  mesoscale  and  laboratory-scale  vortex 
streets  requires  some  theoretical  explanation.   A  well  known  condition  for 
the  formation  of  a  train  of  discrete  vortices  behind  an  obstacle  (usually  a 
circular  cylinder)  in  air  or  water  on  a  laboratory  scale  is  that  the  Reynolds 
number  lie  roughly  in  the  range  102  <  Re  <  105.  Here,  the  Reynolds  number, 
ul/v,    is  defined  in  terms  of  the  ambient  flow  velocity  u,  the  diameter  I   of 
the  obstacle  and  the  kinematic  viscosity  v  of  the  fluid.   If  the  Reynolds 
number  is  lower  than  specified,  the  flow  tends  to  be  laminar;  if  larger,  the 
flow  in  the  wake  of  the  obstacle  tends  to  be  fully  turbulent.   The  periodic 
shedding  of  vortices  represents  the  transition  regime  between  a  laminar  and  a 
turbulent  wake.   On  a  laboratory  scale,  the  required  Reynolds  numbers  are 
consistent  with  the  kinematic  molecular  viscosities  of  the  fluids  used,  but 
if  one  tries  to  explain  the  mesoscale  eddies  in  the  same  way,  kinematic 
viscosities  result  that  are  much  too  large  to  be  attributable  to  molecular 
processes.   Zimmerman  (1969),  Goldman  and  Wilkins  (1973),  and  others  have 
pointed  out  that  the  concept  of  eddy  viscosity  must  be  invoked  if  effective 
Reynolds  numbers  in  the  required  range  are  to  be  obtained.   Eddy  viscosity 
arises  in  a  turbulent  "mixing  region"  between  the  environmental  flow  and 
internal  storm  circulation.  Coefficients  of  eddy  viscosity  of  the  order  of 
102-103  m2  sec-l  are  required  in  the  case  of  thunderstorm-obstacles. 

In  the  theory  of  aeolian  tones  (cf.  Powell,  1964,  and  the  many  refer- 
ences he  cites)  the  vortex-shedding  frequency  (that  is,  the  rate  at  which 
vortices  are  shed  from  one  side  of  the  obstacle)  corresponds  to  the  frequency 
of  the  emitted  tone.   Laboratory  experiments  have  established  empirically 
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that,  within  the  required  range  of  Reynolds  numbers,  the  vortex-shedding 
frequency  fs  is  nearly  proportional  to  the  stream  velocity  U  and  inversely 
proportional  to  the  obstacle  diameter  D  (Figure  11).   The  constant  of  pro- 
portionality is  called  the  Strouhal  number  S,  so  that 


f  =  ~ 


SU 

D 


(6.6) 


Measured  values  of  S  are  always  in  the  neighborhood  of  0.2,  even  when  eddy 
viscosity  plays  a  role  (Goldman  and  Wilkins,  1973).   We  can  thus  estimate  the 
U/D  ratios  that  would  be  required  to  explain  the  observed  infrasound  wave 
periods  t.   Using 


D/U  =  T/5, 


(6.7) 


we  get  for  T  =  20  s,  U/D  =  0.25,  and  for  x  =  200  s,  U/D  =0.025.   If  we  let 
U  be  no  larger  than  50  m  s~l,  then  D  must  be  smaller  than  200  m  (for  t  =  20  s) 
or  2  km  (for  T  =  200  s) .   These  obstacle  sizes  seem  rather  small  in  view  of 
radar  data,  which  suggest  an  obstacle  size  of  5  to  10  km.  A  10  km  diameter 
yields  a  1000  s  shedding  period  (for  u  =  50  m  s~l). 

Using  a  formula  given  by  Phillips  (1956)  for  the  intensity  of  aeolian 
tones,  we  get  for  the  total  radiated  power 


2   2 

n  =  o.2tt  p 


I   DU6S2/C3: 


(6.8) 


where  p  is  the  atmospheric  density  and  I   is  the  obstacle  length.   Taking 
u  =  50  m  s~l 


5  km,  p  =  1.2  kg  m"3,  S  =  0.2,  and  C  =  340  m  s"1, 
we  get  II  =  1.27  x  10y  W.   The  power  level  is  thus  in  the  experimental  range, 
but  the  wave  period  (about  8  min  for  the  values  used)  is  too  long. 


I   =   D  = 
>9 
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Figure  11.      Sketch  of  the  streamlines  in  the  flow  past  a  cylindrical  obstacle 
under  conditions   that  permit  a  train  of  wake  vortices  to  be  shed  from  the 
obstacle.     Acoustic  waves  are  emitted  at  the  shedding  frequency  in  a  dipole- 
like  pattern. 
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7.   HEAT-DRIVEN  OSCILLATIONS 

It  has  been  known  for  a  long  time  that  gas  flows  through  systems  con- 
taining a  heat  source  can  produce  acoustic  oscillations  that  are  coupled  with 
a  periodic  heat  release  by  the  heat  source  (Rijke,  1859;  Rayleigh,  1878, 
1879;  Merk,  1956;  Mawardi,  1956).   Such  oscillations  are  called  heat-driven 
oscillations.  In  general,  coupling  between  a  flow  and  a  heat  source  occurs 
when  fluctuations  in  the  flow  parameters  (pressure  or  temperature,  for  ex- 
ample) influence  the  rate  of  heat  release,  and  when,  in  turn,  fluctuations  in 
heat  release  influence  the  flow  state  variables.   Instability  and  oscilla- 
tions occur  when  the  heat  source  reacts  to  perturbations  in  the  proper  phase 
to  reinforce  the  perturbations.   The  initial  fluid  fluctuations  may  be  either 
small -amplitude  acoustic  waves  or  random  fluctuations  carried  by  the  flow. 
Oscillations  at  particular  frequencies  are  often  enhanced  by  direct  acoustic 
feedback,  as  when  the  system  is  enclosed  in  an  acoustic  resonator.   Such 
oscillations  are  observed  in  certain  kinds  of  "singing  flames,"  and  sometimes 
have  destructive  effects  in  industrial  combustion  systems  and  in  jet  and 
rocket  engine  combustion  chambers  (see  the  review  paper  by  Putnam  and  Dennis, 
1956). 

The  question  I  ask  here  is  whether  similar  oscillations,  and  the  asso- 
ciated sound  radiation,  could  be  caused  by  the  interaction  between  a  convec- 
tive  storm's  updraft  and  the  release  of  the  latent  heat  of  vaporization  where 
water  droplets  condense.   One  could  qualitatively  argue  as  follows: 

Suppose  initial  fluctuations  in  (say)  temperature  pass  through  the 
saturated  air  inside  a  storm,  either  in  the  form  of  small -amplitude  acoustic 
waves  or  as  fluctuations  in  the  updraft  current.   An  air- temperature  increase 
should  quickly  heat  up  the  water  droplets  and  cause  an  increase  in  evapora- 
tion (assuming  buoyancy  effects  occur  much  more  slowly) .  As  the  droplets 
evaporate,  they  cool  and  draw  heat  from  the  surrounding  air,  thus  cooling  it 
by  conduction.   But  this  process  takes  time:   There  are  characteristic  time 
scales  associated  with  the  transfer  of  mass,  heat  and  momentum  among  the 
water  droplets,  the  dry  air  and  the  water  vapor  components  of  the  mixture. 
If  the  air  temperature  responds  to  the  droplet  cooling  just  as  a  negative 
temperature  phase  of  the  initial  fluctuation  passes  through  the  region,  then 
that  frequency  component  of  the  fluctuations  should  be  reinforced.   During 
the  negative  temperature  cycle,  increased  condensation  would  be  encouraged, 
releasing  more  latent  heat,  which  heats  up  the  surrounding  air  in  time  for 
the  arrival  of  the  next  positive  temperature  phase  of  the  initial  fluctua- 
tion.  Thus,  the  frequency  component  of  the  temperature  fluctuations  that 
coincides  with  the  time  scale  of  latent  heat  transfer  would  be  amplifed, 
would  extract  energy  from  the  heat  source,  and  would  radiate  it  away  as 
sound. 

This  qualitative  argument  can  be  submitted  to  quantitative  examination 
by  analyzing  the  propagation  of  plane  acoustic  (or  acoustic-gravity)  waves  in 
an  infinite,  saturated  environment,  as  has  been  done,  for  example,  by  Einaudi 
and  Lalas  (1973)  and  by  the  earlier  workers  they  cite.   Their  results  show 
that,  for  waves  with  very  long,  gravity-wave  periods,  the  effect  of  condensa- 
tion is  essentially  to  change  the  effective  Brunt-Vaisala  frequency  in  the 
wave- dispersion  relation.   At  these  low  frequencies,  the  water  droplet  mass, 
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velocity,  and  temperature  respond  fully  to  environmental  (i.e.,  wave-asso- 
ciated) changes  and  remain  nearly  in  equilibrium  with  the  air.   For  audible 
acoustic  waves,  on  the  other  hand,  the  condensation  time  constants  are  much 
longer  than  a  wave  cycle,  so  that  the  condensation  process  just  barely  gets 
started  before  the  wave  cycle  reverses.   The  net  effect  of  the  resulting  lag 
is  a  relaxation  type  of  absorption  in  saturated  air.   The  frequency  range 
that  interests  us  lies  between  these  two  extremes,  when  a  wave  period  may  be 
comparable  with  the  time  scale  of  the  phase- exchange  process. 

The  problem  of  sound  propagation  in  a  condensing  vapor  has  been  analyzed 
by  Marble  and  Wooten  (1970),  Marble  (1969),  Cole  and  Dobbins  (1970),  Einaudi 
and  Lalas  (1973),  and  others.   The  analysis  by  Marble  and  Wooten  appears  to 
apply  at  any  wave  frequency  and  is  most  readily  applied  to  a  quantitative 
treatment  of  our  problem.   They  have  obtained  wavelike  solutions  of  the  six 
equations  representing  continuity,  motion,  and  energy  for  the  gas  phase, 
motion  and  energy  of  the  liquid  (droplet)  phase,  and  a  vapor  phase  equation 
of  continuity.   They  have  obtained  a  general  acoustic-wave  dispersion  rela- 
tion of  the  form 


CK2    (Z1-iw)   (Z2-iw)   (Z3-ia)) 
oo    {Q ,-ioj)   (STL-ioo)   (ft  -ioo) 


(7.1) 


which  applies  over  all  frequency  ranges.   The  parameters  Z^  and  co^  are  com- 
plicated functions  of  the  vapor  and  droplet  mass  fractions,  the  specific  heat 
ratio  and  the  latent  heat  of  vaporization;  K  =  Kj+il^,  the  complex  wave  num- 
ber.  From  (7.1)  we  obtain  the  imaginary  part  of  the  wave  number,  i.e.,  the 
attenuation  coefficient : 


2CK2   o>(Z  /ft  -1)  •  oo(Z2/^2l)  •  u)(Z  /fi  -1) 

W     (l+oo2/^2)  •  (l+co2/^22)  •  (1+oj2/^32) 


(7.2) 


CK2 

assuming  -  1.   If  amplification  were  possible,  the  attenuation  coeffi- 
cient would  go  negative  for  some  realistic  values  of  the  system  parameters. 
However,  the  quantities  Z^/ft^  are  always  greater  than  unity.   Thus  the  wave 
is  always  attenuated,  and  the  attenuation  coefficient  exhibits  maxima  at  wave 
frequencies  corresponding  to  co=ftj_,  which  are  characteristic  frequencies 
related  to  the  mass,  momentum,  and  heat- transfer  times  of  the  system. 

For  wave  frequencies  much  smaller  than  the  frequencies  corresponding  to 
the  usual  thermal  and  viscous  time  constants  in  the  system  (i.e.,  about  100 
Hz)  Marble  and  Wooten  get  an  attenuation  coefficient  a 

-  -  j   K  (n-1)   j     ,  (7.3) 

1  +  (wxD/Kp) 
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where  co  is  the  acoustic  wave  frequency,  y   is  the  specific  heat  ratio  and  IC^ 
and  K„  are  the  mass  fractions  of  water  vapor  and  water  droplets  (assumed 

small);  n.  is  the  latent  heat  parameter  L/CpT,  where  L  is  the  latent  heat  of 
vaporization,  and  Cp  is  the  specific  heat  at  constant  pressure  of  the  gas  and 
T  is  its  temperature;  Tq  is  the  mass- diffusion  time  associated  with  droplet 
evaporation.   Thus,  attenuation  at  these  frequencies  always  increases  with 
the  vapor  mass  fraction  and  the  latent  heat  parameter.   Attenuation  reaches  a 
maximum  at  a  frequency  co  =  Kp/xn,  which,  for  typical  cloud  parameters,  is 
about  a  1-s  wave  period. 

To  reconcile  this  result  with  the  qualitative  argument  at  the  beginning 
of  this  section,  it  is  useful  to  examine  the  transient  solutions  obtained  by 
Marble  and  Wooten,  which  reveal  the  response  of  the  system  variables  to  a 
step  discontinuity  in  air  pressure  and  temperature.   The  water  droplets  res- 
pond in  two  phases  to  an  initial  step  increase  in  air  temperature:   First  the 
droplets  are  exponentially  warmed  by  conduction  from  the  air  on  a  time  scale 
of  about  0.01  s,  with  relatively  little  change  in  air  temperature.   Evapora- 
tion begins  when  the  saturation  vapor  pressure  near  each  droplet  surface  has 
been  raised  above  the  ambient  vapor  pressure.   It  takes  about  0.01  s  for 
evaporation  to  begin  after  an  initial  air  temperature  rise.   Evaporation 
cools  the  droplets  and  the  newly  evaporated  vapor,  and  the  gas  is  eventually 
cooled  exponentially  by  conduction  on  a  time  scale  of  about  1  s.   The  final 
equilibrium  gas  temperature  is  determined  by  the  amount  of  liquid  present  and 
the  latent  heat  of  vaporization.   The  net  effect  is  that  droplet  evaporation 
absorbs  heat  from  the  gas  whenever  the  temperature  perturbation  is  positive; 
conversely,  condensation  at  the  droplet  surface  releases  heat  whenever  the 
temperature  perturbation  is  negative.   Whatever  the  sign  of  the  air  tempera- 
ture perturbation,  the  heat  source  acts  to  oppose  it.   Thus,  the  heat  source 
responds  to  a  sinusoidal  temperature  perturbation  whose  frequency  corresponds 
to  the  evaporation  time  scale  (about  1  s) ,  by  opposing  the  excitation  with  a 
slight  time  lag.   The  net  result  is  absorption  of  the  perturbing  wave,  and 
not  reinforcement.   Thus,  the  mechanism  is  not  a  viable  source  process. 


7.1  Condensation  in  a  Supersaturated  Environment 

Another  acoustic  source  mechanism  has  been  suggested  that  is  related  to 
the  heat-driven  oscillations  just  examined.   If  condensation  commences  in  a 
supersaturated  region  of  the  atmosphere,  the  condensation  boundary  should 
move  outward  with  a  speed  limited  by  the  sound  speed.   Corresponding  density 
changes  should  move  with  the  boundary,  and  an  acoustic  impulse  should  be 
radiated.   Rapidly  building  towers  in  some  cumulonimbi  may  be  visual  manifes- 
tations of  such  moving  condensation  boundaries. 

It  is  difficult  to  construct  a  realistic  quantitative  model  of  this 
process  from  which  the  acoustic  radiation  could  be  estimated.   The  process  is 
mentioned  here  only  for  completeness.   However,  the  radiation  should  qualita- 
tively resemble  the  impulsive  emissions  calculated  for  the  imploding  sphere 
in  Section  5.2  on  electrostatic  sound.   Similar  objections  could  be  raised 
about  the  impulsive  nature  of  such  a  source. 
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8.   SUMMARY  AND  CONCLUSIONS 


8. 1  Relative  Merits  of  the  Models 

The  simple-source  models  predict  the  radiation  of  a  pressure  pulse  of 
about  the  same  duration  as  the  growth  phase  of  a  convection  cell  --  quite 
different  in  character  from  the  observed  infrasound.   There  is  no  experi- 
mental evidence  that  fluctuations  on  shorter  time  scales  emit  coherently  over 
the  spatial  extent  of  a  storm.   However,  fluctuations  in  the  amplitude  of 
infrasound  measured  both  in  the  ionosphere  and  on  the  ground  do  exhibit  a 
tendency  toward  a  20  to  30  min  periodicity  which  has  previously  been  attri- 
buted to  storm  life  cycle  (Georges,  1973). 

Fluctuations  on  smaller  spatial  scales  would  emit  sound,  but  their 
contributions  would  add  incoherently  over  the  extent  of  the  storm.   The 
strength  and  spectral  content  of  such  emissions  are  best  modeled  stochas- 
tically: A  turbulent- emission  model  predicts  an  average  radiated  power  spec- 
tral density  proportional  to  f~7/2}  with  most  of  the  emitted  power  concen- 
trated at  frequencies  corresponding  to  the  outer  scale  of  the  turbulence. 
Realistic  velocity  fluctuations  of  25  to  40  m  s~l  extending  over  a  10  km- cube 
storm  volume  would  radiate  a  total  broadband  power  of  10^  to  10^  W.   The 
spatial  and  temporal  scales  of  these  fluctuations  (the  turbulent  outer  scale) 
would  cause  the  radiation  to  be  strongest  at  20  to  30  min  wave  periods  and 
considerably  weaker  at  the  tens  of  seconds  periods  observed.   Also,  such  a 
model  cannot  explain  the  observed  trains  of  narrow-band  infrasound  sometimes 
observed. 

The  thermo-mechanical  mechanisms  appear  not  to  be  applicable  to  severe 
storms,  because  the  heat-release  mechanism  (condensation)  is  always  out  of 
phase  with  wave- associated  temperature  fluctuations. 

The  electrical  models  appear  to  provide  insufficient  power  in  the  ob- 
served spectral  region;  furthermore,  it  is  difficult  to  explain  the  nearly 
continuous  wave  trains  observed  in  terms  of  an  essentially  impulsive  source. 
Enough  uncertainty  remains,  however,  that  an  experimental  test  would  be 
desirable. 


8.2  A  "Most  Likely"  Model 

The  close  association  of  infrasound  with  tornadic  storms  and  the  degree 
with  which  the  theoretical  model  of  vortex  radiation  matches  the  observations 
make  the  vortex-sound  mechanism  the  "most  likely"  candidate.  The  observation 
that  the  infrasound  seems  to  be  related  to  tornado  formation,  but  in  some 
cases  precedes  it  (Georges  and  Greene,  1975),  along  with  the  calculations  of 
the  vortex  size  required  to  match  the  observed  radiated  power  and  frequency, 
suggest  that  the  rotating  columns  of  air  inside  some  thunderstorms  called 
"tornado  cyclones"  may  be  the  infrasound  emitters.   Radar  observations  show 
that  tornado  cyclones  form  up  inside  storms  as  much  as  an  hour  before  tor- 
nadoes touch  ground.   (Note  that  Georges  and  Greene  (1975)  report  that 
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infrasonic  emissions  often  precede  observed  tornadoes  by  up  to  an  hour.) 
Evidence  from  the  deviation  and  curvature  of  tornado  tracks  (especially 
multiple  tornadoes)  suggest  that  the  tornadoes  may  be  satellite  vortices  of 
the  moving  longer  tornado  cyclone,  as  suggested  by  Agee  et_  al.  (1976).   But 
whether  or  not  tornadoes  that  reach  the  ground  actually  form  from  a  given 
meso-cyclone,  the  infrasound  radiating  mechanism  should  operate. 

However,  the  quasi-random  character  that  some  of  the  observed  emissions 
assume  suggests  that  some  features  of  the  turbulent  source  model  are  also 
relevant.   The  two  models  are  not  incompatible.   A  "hybrid"  model,  combining 
features  of  each,  in  fact  appears  to  reproduce  the  most  features  of  the 
observations.  The  actual  radiation  process  could  involve  up  to  several  indi- 
vidual vortex  radiators  at  various  stages  of  growth  and  decay  inside  a  single 
storm  complex  (but  not  necessarily  in  a  single  storm  cell).   Many  times,  as 
is  observed,  a  single  source  could  dominate,  radiating  nearly  monochromat- 
ically.   At  other  times,  the  emissions  may  be  the  superposition  of  many 
sources.   The  key  hypothesis  that  joins  the  two  mechanisms  is  that  the  two 
are  fundamentally  the  same  physical  process,  but  are  only  described  differ- 
ently.  A  very  large  number  of  randomly  distributed  eddies  or  vortices  ap- 
proaches what  would  normally  be  called  turbulence,  and  their  total  radiation 
could  only  be  computed  statistically.   Averaged  over  many  storm  events,  then, 
the  radiation  would  obey  the  spectral  predictions  of  the  turbulent  model;  for 
example,  the  average  radiated  power  spectrum  would  follow  an  f-7/2  law  down 
to  wave  periods  corresponding  to  storm-scale  motions,  i.e.,  20  to  30  min.   At 
the  other  extreme,  a  single  vortex  that  degenerates  in  an  orderly  way  can  be 
modeled  analytically  and  its  radiation  calculated  exactly.   The  situation 
that  we  encounter  in  nature  seems  to  lie  between  these  two  extremes,  where 
neither  model  is  completely  satisfactory,  but  in  different  senses,  each  is. 
This  aspect  of  the  model  needs  to  be  refined  by  being  placed  on  firm  analy- 
tical ground. 

Application  of  the  simple  model  for  the  acoustic  radiation  from  a  coro- 
tating  vortex  pair  to  the  problem  of  thunderstorm  radiation  provides  a  bridge 
between  two  different  approaches  to  the  problem.   On  one  hand,  there  are 
completely  deterministic  models  of  coherent  acoustic  emissions  based  on 
certain  storm  processes  known  to  resemble  fluid  monopole  and  dipole  sources; 
on  the  other  hand,  it  is  possible  to  model  the  flow  inside  a  storm  stochas- 
tically, that  is,  as  "turbulence"  and  to  estimate  the  acoustic  radiation 
according  to  Lighthill's  model  based  on  a  random  distribution  of  acoustic 
quadrupoles.   The  disadvantage  of  the  former  kind  of  model  is  that  most  of 
the  phenomena  inside  storms  that  we  know  about  and  can  model  determinis- 
tically  lead  to  the  radiation  of  impulses  only,  or  at  least  lead  to  very  much 
longer  period  radiation  than  we  require.   The  disadvantage  of  a  stochastic 
approach  is  that  it  masks  most  of  the  physics  that  would  relate  storm  pro- 
cesses to  the  emissions.   The  vortex-sound  model  provides  what  is  needed:   a 
basic  radiating  mechanism  that  can  predict  either  coherent  or  incoherent 
emissions  (both  are  observed) ,  depending  on  the  size  and  number  of  individual 
emitters  inside  a  storm  and  their  lifetimes. 


49 


8.3  Some  Observational  Tests 

The  electrical  mechanisms  lend  themselves  to  tests  that  observe  simulta- 
neously the  infrasonic  and  electromagnetic  energy  that  storms  radiate.  For 
example,  if  the  infrasound  is  part  of  the  acoustic  remnant  of  the  cylindrical 
shock  generated  by  lightning  strokes,  then  there  should  be  a  spatial  and 
temporal  correlation  between  the  infrasound  sources  and  the  sources  of  "atmo- 
spherics," the  radio- frequency  electromagnetic  energy  radiated  by  lightning. 
Such  a  test  has  been  performed  and  is  reported  by  Beasley  et  al.  (1976),  with 
negative  results.   Their  results  also  show  that  other  electromagnetic  mech- 
anisms (like  the  electrostatic-relaxation  process)  are  also  unlikely  explana- 
tions for  the  observed  infrasound. 

The  vortex- sound  mechanism  suggests  several  tests.   Again,  simultaneous 
mapping  of  an  emitting  storm's  velocity  field  with  high-resolution  Doppler 
radar  would  permit  a  direct  association  between  the  radiation  and  the  dyna- 
mics and  structure  of  mesoscale  vortices  inside  storms.   The  several  Doppler 
and  dual-Doppler  meteorological  radars  now  coming  into  operation  increase  the 
chances  that  such  simultaneous  observations  will  eventually  happen  (Figure  12) 
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Figure  12.     A  map  of  horizontal  air-veto oity  vectors  at  an  altitude  of  6  km 
above  ground  inside  a  Colorado  thunderstorm.      The  velocity  field  was  mapped 
with  the  WPL  dual-Doppler  radar   (Kropfli  and  Miller 3    1976).      The   length 
(and  intensity)   of  each  vector  is  proportional  to  aiv  speed,   with  mean  re- 
moved,   at  that  location,   with  the  highest  speeds  reaching  10  m  s~l.     A 
strong  vortex  is  visible  at  lower  left,   and  two  weaker  ones  at  the  top  and 
at  the   lower  right  can  also  be  seen.      The  radar  resolution  does  not  allow 
the  sub- vortex  structure  to  be  seen  with  any  detail. 
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If  the  infrasound  is  radiated  from  vertical  vortices,  the  radiation 
should  have  directional  characteristics  that  could  be  observed  with  sensors 
spaced  azimuthally  around  the  source.   Azimuthal  (rotational)  symmetry  in  the 
source  would  be  reflected  in  radiation  that  is  nondirectional,  except  possi- 
bly for  small  time  delays  that  may  be  detectable.   Variability  of  propagation 
time  would  make  this  difficult,  however.   Azimuthal  vortex  asymmetries,  such 
as  those  observed  by  Agee  et  al.  (1975)  could  cause  azimuthal  variations  in 
radiated  frequency  that  may  follow  a  predictable  (and  observable)  pattern. 
Frequency- time  spectra  of  emissions  recorded  at  many  observatories  spaced  in 
azimuth  could  reveal  such  patterns. 

Another  test  has  been  devised  to  measure  the  radiation  from  a  more 
accessible  form  of  atmospheric  vortex:   dust  devils.   There  is  photographic 
evidence  that  dust  devils  exhibit  multiple-vortex  and  turbulent  structures 
(Figure  13).  Calculations  using  the  formulas  of  Section  6.3.1  show  that  dust- 
devil  emissions  should  be  detectable  with  existing  infrasonic  microphones. 
We  attempted  such  an  experiment  but  found  no  emissions  that  could  be  attri- 
buted unambiguously  to  several  large  dust  devils  that  passed  near  our  micro- 
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Figure  13.  Photograph  (by  the  author)  of  a  dust-devil  near  Phoenix,  Ariz.3 
showing  helical  structure.  Motion  pictures  of  the  same  dust-devil  reveal 
satellite  vortices  embedded  in  the  flow. 
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phones.  However,  man-made  interference,  low  equipment  sensitivity  and  the 
small  number  of  dust  devils  sampled  could  easily  explain  our  negative  result 
More  extensive  experiments  should  be  done. 

If  monopole  or  dipole  oscillations  of  storm  scale  existed  with  suffi- 
cient intensity  and  at  the  right  frequency  to  explain  the  infrasound  obser- 
vations, we  would  expect  to  observe  corresponding  large  near- field  (non- 
radiated)  infrasonic  pressure  fluctuations  in  the  vicinity  of  the  emitting 
storms.  It  should  be  straightforward  to  estimate  the  expected  near- field 
intensity  and  then  either  to  examine  existing  pressure  records  (most  of 
which,  however,  may  not  respond  at  the  infrasonic  frequencies)  or  to  devise  < 
simple  sensor  to  detect  them. 
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